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Abstract 

Researchers  att  he  A  irF  orce  Institute  o  fT  echnology  (AFIT)  an  d  the  O  perationally 
Responsive  S  pace  (  ORS)  O  ffice  ha  ve  c  onducted  e  xtensive  vi  bration  testing  a  nd  s  tructural 
modeling  on  the  first  ORS  Plug-and-Play  Satellite  (PnPSAT  I).  The  intent  of  this  research  effort 
is  t  o  ev  aluate  t  he  p  remise  t  hat  cu  rrent  p  ost-integration  s  pacecraft  en  vironmental  t  est 
requirements  can  be  reduced  or  modified  using  accurately  tuned  finite  element  (FE)  models.  As 
part  of  this  research,  modal  testing  was  conducted  on  the  PnPSAT  I  structural  panels  at  AFIT. 
The  modal  te  sting  was  part  o  f  a  m  uch  1  arger  s  cries  o  f  ex  perimental  t  rials  o  n  v  arious 
configurations  of  PnPSAT  I  at  the  Air  Force  Research  Laboratory  (AFRL)  facilities  at  Kirtland 
Air  Force  Base  (KAFB).  Multiple  sets  of  vibration  data  were  also  collected  from  accelerometers 
on  PnPSAT  I  from  standard  and  modified  spacecraft  prelaunch  sine  sweep  and  random  vibration 
tests.  The  modal  data  collected  at  AFIT  is  used  to  tune  two  PnPSAT  I  panel  FE  models  and  the 
random  vi  bration  da  ta  collected  at  K  AFB  is  used  to  tu  ne  th  e  c  omplete  s  atellite  f  or  o  ne 
configuration.  T  he  goal  i  s  t  o  cr  eate  an  a  ccurate  F  E  m  odel  c  apable  o  f  p  redicting  t  he  d  ynamic 
response  in  a  f  requency  range  of  0-300  Hz  of  various  PnPSAT  configurations.  T  his  modeling 
and  tuning  effort  will  be  validated  by  comparing  FE  model  predictions  with  measured  vibrational 
response  from  the  previously  mentioned  experimental  trial. 
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ACCURATE  DYNAMIC  RESPONSE  PREDICTIONS  OF  PNPSAT  I 


I.  Introduction 

Satellites  are  very  complex  systems  that  must  meet  intensive  and  strict  design,  manufacture, 
test,  and  integration  requirements  for  a  successful  mission.  Rigorous  post-integration  acceptance 
tests  a  re  c  ritical  to  guard  a  gainst  failures  during  1  aunch  and  de  ployment.  Presently,  pos  t- 
integration  satellite  environmental  tests,  such  as  vibration  and  thermal  vacuum  (TVac),  typically 
exceed  several  months  in  duration  for  traditional  satellites  Wertz  [1999],  Yee  [2005],  However, 
Plug-and-Play  (PnP)  satellites  are  designed  to  fulfill  rapid,  time  sensitive  mission  requirements. 
Because  satellites  are  typically  very  expensive,  traditional  satellites  are  mostly  designed  as  a  one- 
of-a-kind  f  light  a  rticles  and  i  ts  components  a  re  or  dered  or  bui  It  i  n  as-required  qua  ntities. 
Confidence  i  n  t  he  f  ully-integrated  s  atellite  i  s  onl  y  gained  t  hrough  e  xhaustive  pr  e-flight 
acceptance  t  ests.  C  onversely,  P  nP  s  atellites  ar  e  co  nstructed  f  rom  a  wide  ar  ray  o  f  k  nown 
components  and  are  designed  to  minimize  integration  issues.  L  ike  components  on  t  raditional 
satellites,  all  PnP  components  will  have  passed  individual  and  possibly  subsystem  ope  rational 
acceptance  te  sts.  P  rior  to  s  atellite  c  onstruction,  after  the  components  h  ave  be  en  s  elected  f  or 
flight,  additional  and  possibly  exhaustive  environmental  tests,  including  hardware  and  software, 
are  typically  conducted  to  significantly  increase  confidence  in  the  fully-integrated  satellite. 

These  en  vironmental  t  ests  ar  e  n  ecessary  t  o  en  sure  t  hat  t  he  s  atellite  i  s  cer  tified  flight-ready 
and  that  no  issues  or  anomalies  are  overlooked.  In  general,  all  flight  hardware  is  exhaustively 
tested  i  n  or  der  toe  ertify  t  hat  i  t  c  an  s  urvive  1  aunch  a  nd  ope  rate  i  n  t  he  s  pace  e  nvironment.  In 
order  to  reduce  risk,  full  assembly  integration  and  test  (AIT)  program  is  conducted  on  each  level 
of  a  ssembly  ranging  f  rom  c  omponents  a  nd  s  ubsystems  t  o  t  he  e  ntire  s  atellite.  However, 
traditional  timelines  to  complete  post-integration  vibration  and  TVac  tests  must  be  reduced  for 
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PnP  s  atellites  t  o  be  1  aunched  i  n  da  ys  a  s  oppos  ed  t  o  m  onths.  I  n  or  der  t  o  de  termine  i  f  t  est 
processes  c  an  be  r  educed  or  m  odified  w  bile  m  aintaining  hi  gh  reliability,  A  FIT  and  O  RS 
researchers  are  conducting  experiments  to  evaluate  traditional  acceptance  tests,  reduced,  and  new 
acceptance  tests  on  a  fully- integrated  satellite,  this  set  of  experiments  is  also  referred  to  as  the 
Rapid  AI&T  Demonstration. 


Figure  1 :  PnPSAT  fully  configured  ready  for  test. 

1.1  ORS  and  PnP  Test  Program 

The  ORS  office  is  responsible  for  the  development  of  capabilities  that  provide  assured  space 
power  focusing  on  t  imely  mission  c  onstraints  i  n  or  der  t  o  s  atisfy  t  he  J  oint  F  orce  C  ommanders  ’ 
needs  by  employing,  deploying,  or  developing  new  systems  to  sustain  or  augment  capabilities 
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Bhopale  [2009].  ORS  consists  of  three  distinct  tiers  eaeh  responsible  for  a  different  task,  these 
tiers  are  shown  in  Figure  2. 


Tier-1 


Tier-2 


Tier-3 


-  Rapidly  exploit  existing  capabilities 
•Net-Cenb'ic 

■  Open  Architecture 

-Data 

-  Disccverabte 

-  Subscribeble 
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-  Interoperable 


Replenish 
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ORS  Office 


ORS  Activities 

Figure  2:  The  three  tiers  of  the  ORS  Office  Bhopale  [2009], 


Tier-1  is  comprised  of  capabilities  readily  available  that  provide  highly  responsive  space 
effects.  Therefore,  these  solutions  utilize  existing  space  systems,  ground  systems,  operations  and 
processes  an  d  t  ypically  are  em  ployed  w  ithin  a  f  ew  h  ours-to-days.  If  T  ier- 1  i  s  una chievable  a 
Tier-2  solution  is  considered.  T ier-2  solutions  utilize  field-ready  capabilities  which  exist  in  the 
days-to-weeks  time  frame.  Consequently  Tier-3  exists  to  create  technologies  that  will  eventually 
be  classified  as  Tier-2  once  they  are  field  ready.  For  a  more  detailed  description  of  ORS  refer  to 
Bhopale  [2009].  P  nPSAT  I  falls  into  the  Tier-2  classification,  PnPSAT  I  is  meant  to  be  fully 
assembled  and  ready  for  mission  within  6  days  of  being  called  upon. 

Since  the  beginning  of  satellite  assimilation  and  deployment,  satellite  design  has  been  an 
extremely  non- standardized  and  complex  process.  U  sually,  hi ghly  experienced  personnel  are 
required  to  synthesize,  create,  build,  test,  and  launch  a  satellite.  Like  the  advent  of  the  personal 
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computer  (P  C)  with  the  invention  of  uni  versal  serial  bus  ( USB),  strides  ha  ve  be  en  t  aken  t  o 
simplify  th  e  way  satellites  ar  e  co  nstructed.  T  he  A  ir  Force  R  esearch  Lab,  S  pace  Vehicles 
Directorate  (ARFL/RV)  has  undertaken  two  efforts  in  small  satellite  development  to  demonstrate 
the  effectiveness  of  the  technology  and  show  it  is  beneficial.  The  first  effort,  PnPSAT  1  utilizes 
the  space  Plug-n-Play  avionics  (SPA)  interface  standard.  Implementation  of  this  SPA  interface 
standard  has  proven  that  rapid  development,  integration  and  test  is  possible.  In  the  second  effort, 
PnPSAT  I  uses  the  next  generation  of  SPA  components  for  a  larger  bus  focused  on  ORS  needs  to 
ensure  custom  performance  at  commodity  prices.  In  short,  all  the  components  exist  and  are 
meant  to  be  like  the  name  suggests,  Plug-n-Play,  similar  to  the  PnP  standard  for  PCs.  The  idea  is 
that  t  he  us  er  can  pi  ug  c  omponents  i  nto  t  he  satellite  us  ing  s  tandard  c  onnectors  across  a  11 
components  and  quickly  assemble,  test,  and  launch  a  satellite  in  order  to  replenish  or  populate  a 
constellation  to  conduct  a  mission.  The  predicted  timeline  for  ORS  Tier-2  satellites  allows  only 
a  three  day  window  for  integration  of  testing  (Rapid  AIT).  Consequently,  leaving  only  two  days 
to  fully  integrate  and  launch  the  satellite. 

As  p  art  o  f  t  he  a  forementioned  R  apid  A  IT,  it  is  imp  ortant  to  minimize  vibe  a  nd  t  hermal 
testing.  T  hese  expensive  t  ests  a  re  m  anpower  a  nd  time  consuming;  elimination  o  r  a  t  le  ast  a 
reduction  of  these  tests  is  of  critical  importance  to  the  ORS  mission.  T  he  overall  test  program 
for  reducing  these  post-integration  is  described  below. 

1.2  Overall  Rapid  AI&  T  Program 

The  Rapid  AIT  demonstration  is  designed  to  explore  the  premise  that  current  spacecraft  test 
requirements  can  be  reduced  or  modified  for  ORS  PnP  satellites  by  modifying  processes,  testing 
strategies,  and  personnel  dependencies  Baghal  [2010].  There  are  four  trials  in  the  Rapid  AI&T 
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Demonstration  pi  an.  E  ach  trial  is  designed  to  investigate  one  variables’  dependency  on  that 
same  variable.  F  or  example,  in  Trials  1  &  2,  the  only  difference  is  the  personnel  performing  the 
assembly  and  test.  The  two  groups  were  referred  to  as  Team  A  and  Team  B.  Team  A  consisted 
of  personnel  intimately  familiar  w  ith  P  nPS  AT.  T  earn  B  consisted  of  spacecraft  t  echnologists 
that  are  not  familiar  with  PnPSAT  specifically.  Team  B  trained  during  Trial  1  and  performed  the 
assembly  and  test  on  T  rial  2.  T  he  updates  to  AI&T  flow  and/or  s pacecraft  c onfiguration  for 
Trials  3  and  4  will  be  determined  upon  an  assessment  from  Trials  1  and  2. 

After  each  PnP  satellite  disassembly  and  reassembly,  traditional  acceptance  tests,  reduced, 
and  new  acceptance  tests  are  conducted  on  the  components  and/or  subsystems.  Dynamic  and 
thermal  m  odels  ar  e  co  rrelated  w  ith  t  he  co  llected  d  ata  an  d  p  redictions  will  b  e  m  ade  o  n  the 
response  of  the  fully-integrated  satellite.  Later  tests  include  various  configurations  in  order  to 
determine  to  the  capability  of  models  to  predict  responses  of  various  configurations.  A  fter  the 
data  is  analyzed,  th  e  r  esults  ma  y  in  dicate  th  at  p  ost-integration  t  ests  can  b  e  r  educed  w  hile 
incurring  little  increase  in  risk. 

1.3  Finite  Element  Modeling 

In  order  to  further  reduce  rigorous  post-integration  testing,  a  finite  element  (FE)  model  of 
PnPSAT  I  was  created  by  Spaceworks.  T  he  FE  model  was  created  from  a  CAD  model  from 
which  the  satellite  was  made  from.  This  CAD  model  is  extremely  accurate  because  it  includes 
even  the  finer  details  such  as  every  bolt,  washer  and  screw,  as  seen  in  Figure  3.  From  this  CAD 
model,  t  he  m  ass  mo  ments  o  f  in  ertia  (MOI)  of  the  s  atellite  can  b  e  es  timated  which  ar  e  u  sed 
primarily  for  launch  vehicle  integration. 
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Figure  3:  Fully  deployed  high  fidelity  model  of  PnPSAT  1. 


Spaceworks,  the  company  that  created  the  original  CAD  model  was  tasked  to  create  an  FE 
model  of  PnPSAT  I.  T  heir  FE  model  was  used  to  predict  the  stmcture’s  natural  frequencies. 
The  FE  model  has  368,850  nodes  and  246,295  elements  and  1,104,633  DOF,  see  Figure  4. 


Figure  4:  (a)  Spaceworks  high  fidelity  PnPSAT  I  model,  (b)  close-up  view  in  top  of  Z  panel 
showing  rigid  links  (dark  thick  lines),  represented  by  point  masses  with  links,  these  links  connect 
to  the  point  masses  representing  the  components  to  the  satellite's  panel. 
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The  Spaceworks  FE  model,  although  very  refined,  has  too  many  FE  DOF  for  tuning  in  a 
reasonable  amount  of  time.  Tuning  involves  minimizing  the  differenee  between  analytieal  and 
measured  modes  while  eonstraining  the  natural  frequencies.  In  order  to  tune  PnPSAT  1,  a  new 
FE  m  odel  w  ith  s  ignificantly  fewer  el  ements  n  ceded  t  o  b  e  cr  eated.  Like  the  Spaceworks  F  E 
model,  the  Spaceworks  CAD  model  served  as  the  starting  point  for  the  new  reduced  FE  model, 
development  and  tuning  is  the  focus  of  this  research. 

1.4  Research  Overview 

The  PnP  concept  on  which  PnPSAT  I  was  initially  conceived  was  aimed  at  the  elimination  or 
reduction  of  e  xpensive  a nd  t  ime  c  onsuming  A  IT.  T  he  O  RS  m  ission  i  s  t  o  1  aunch  pa  yloads  t  o 
space  in  less  than  six  days.  T  he  objective  of  this  research  is  to  investigate  th  e  p  ossibility  o  f 
creating  accurate  F  E  m  odels  tom  odel  P  nPS  AT  lino  rder  top  redict  i  ts  n  atural  f  requencies  i  n 
different  configurations.  This  would  eliminate  the  amount  of  vibration  testing  required  making  it 
possible  to  launch  in  six  days. 

The  process  is  initialized  by  generating  FE  models  of  both  PnPSAT  Ts  -Y  and  -Z  panels  (only 
two  panels  were  required  for  modal  testing).  The  FE  models  possess  the  very  similar  geometry 
and  mass  as  the  actual  panels  of  the  satellite.  N  ext,  modal  data  is  extracted  from  the  impulse 
response  from  both  the  -Y  and  -Z  panels.  T  he  FE  models  are  then  tuned  using  the  collected 
modal  data  by  only  changing  the  panels  modulus  of  elasticity.  Asa  result,  the  tuned  FE  model 
panels  are  copied  and  combined  to  create  the  main  structure  of  the  PnPSAT  I.  P  nPSAT  I  FE 
model  structure  is  rigidly  attached  via  coincident  nodes,  the  components  are  attached  to  match 
the  first  test  configuration  in  the  Rapid  AIT  Demonstration.  This  complete  satellite  FE  model  is 
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tuned  us  ing  onl  y  na  tural  f  requency  da  ta  collected  a  t  one  of  t  he  pr  eviously  m  entioned  K  AFB 
shaker  tests. 


A  flow  chart  of  this  tuning  process  implemented  in  this  research  is  in  Figure  5. 


Figure  5:  Overall  process  for  creating  and  tuning  PnPSAT  I FE  model. 
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1.5  Thesis  Organization 


The  second  of  five  chapters  provides  a  literature  review  of  previous  efforts  to  tune  satellite  FE 
models,  t  echnical  ba  ckground  be  hind  F  E  m  odeling/analysis,  c  ollecting  m  odal  da  ta,  a  nd  t  he 
optimization  process.  The  third  chapter  describes  the  methods  used  to  create  the  FE  model,  the 
collection  of  modal  data,  and  tuning  the  models.  Chapter  Four  provides  the  results  of  the  modal 
analysis  t  o  i  nclude  unt  uned,  t  uned  m  odels  a  nd  the  c  omparison  of  results  be  tween  di  fferent 
configurations  a  nd  e  xperimental  da  ta.  T  he  f  inal  c  hapter  c  oncludes  w  ith  a  di  scussion  on  t  he 
conclusions  drawn  from  the  results  and  recommendations  for  any  future  work. 
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II.  Background 


The  f  ollowing  chapter  provides  de  tails  i  nto  t  he  F  E  m  odel.  T  he  f  irst  s  ection  provides  a 
literature  review  of  FE  modeling  and  tuning  processes.  The  second  section  describes  FE  analysis 
to  i  nclude  h  ow  t  o  generate  a  FE  m  odel  and  t  he  a  ssumptions  us  ed.  T  bird,  eigenanalysis  is 
described  starting  from  the  equations  of  motion  (EOM)  through  the  eigenvalue  problem  (EVP) 
ultimately  1  eading  t  o  t  he  m  odal  f  requency  r  espouse  pr  oblem.  The  ne  xt  s  ection  d  etails  th  e 
scanning  laser  vibrometer  used  to  collect  the  panel  impulse  responses  at  hundreds  of  points  to 
include  a  discussion  on  how  the  software  works.  Lastly,  the  overall  process  for  tuning  the  FE 
model  with  the  acquired  modes  and  natural  frequencies  is  discussed. 

2.1  Eigenanalysis  and  FE  Model  Tuning  References 


Eigenanalysis  a  nd  FE  model  t  uning  ha  s  be  en  a  ne  cessity  i  n  t  he  space  1  aunch  a  rena  f  or 
sometime.  An  innovative  approach  was  taken  by  Cobb  [1996]  in  his  attempt  to  model  and  tune  a 
6-m  flexible  frame  structure  from  0  to  70Hz.  A  Ithough  a  different  structure,  a  similar  process 
was  used  to  tune  this  model.  The  same  objective  function  as  presented  later. 


j=Za.  1-i  1-4 


was  utilized,  a  long  w  ith  th  e  s  ame  mo  de  tr  acking  p  rocess.  D  ifferent  s  oftware  was  u  tilized  to 
implement  the  process,  ASTROS  (Automated  Structural  Optimization  System  Software)  versus 
NASTRAN.  T  he  s  tructure’s  m  easured  d  ata  w  as  t  aken  w  ith  acc  elerometers  r  eading  o  nly  o  ne 
direction  and  the  input  excitation  to  the  truss  was  two  linear  momentum  exchange  actuators. 
Mounting  plates  were  developed  to  attach  the  excitation  devices  atop  two  longerons  on  the  free 
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end  of  the  truss.  Driving  the  actuators  in  phase  allowed  for  measurement  of  the  bending  modes 
in  one  axis.  The  torsional  modes  were  excited  by  driving  the  two  180  deg  out  of  phase. 

The  FE  model  was  comprised  of  four  different  types  of  beam  elements  for  a  total  of  100 
elements.  S  imilar  to  the  components  in  PnPSAT  I  FE  model,  the  actuators  and  mounting  plate 
were  modeled  as  lump  masses  and  therefore,  included  in  the  dynamics.  The  vibration  test  was 
conducted  us  ing  0-100  Hz  r  andom  noi  se  s  equence.  F  requency  a  veraged  t  ransfer  f  unctions 
between  the  input  excitation  and  the  eight  accelerometers  were  measured.  The  inverse  discrete 
Fourier  transform  of  t  hese  yielded  t  he  i  mpulse  r  esponse  f  unctions  w  hich  w  ere  i  nput  into  a  n 
eigensystem  realization  algorithm  (ERA).  ERA  is  based  on  extracting  a  state-space  model  from 
the  impulse  response  of  a  system.  For  a  more  detailed  look  at  this  method  refer  to  Juang  and 
Pappa  [1985].  ERA  extracted  modes  and  natural  frequencies  from  the  experimental  data. 

Using  this  measured  data  the  FE  model  was  tuned  using  ASTROS  and  various  F ortran  77 
subroutines  created  f  or  e  igenvector  nor  malization,  m  ode  s  witch  t  esting,  ei  genvalue  an  d 
eigenvector  s  ensitivities,  obj  ective  function  e  valuation  a  nd  pr  ocessing  d  ata  t  o  a  nd  f  rom  t  he 
optimization  modules.  The  final  results  of  the  method  showed  the  potential  of  this  process  and  it 
was  a  good  baseline  for  any  FE  model  tuning.  Overall,  the  same  tuning  process  was  utilized,  and 
the  differences  arose  between  the  methods  of  collecting  and  processing  the  data,  and  the  software 
used  to  tune. 

A  unique  approach  was  used  by  Kammer  [1999]  in  studying  the  structural  integrity  of  Russia's 
Mir  space  station.  M  ir  provided  a  rare  opportunity  to  evaluate  methods  of  detecting  structural 
damage  on  a  structure  currently  on  or  bit.  M  odal  identification  was  used  to  gather  data  on  t  he 
potential  damage  on  Mir;  however,  it  provided  a  formidable  challenge  for  several  reasons.  Mir 
modal  data  was  taken  during  short  time  duration  docking  events;  also,  it  was  coupled  with  the 
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availability  of  only  a  few  data  sets  whieh  exeluded  the  effeetiveness  of  averaging.  In  addition, 
due  to  the  docking  element  it  makes  the  problem  nonlinear  making  the  use  of  ERA  unsuitable  for 
this  problem.  T  ypically,  during  a  modal  test,  the  input  and  output  are  measured  from  a  pulse, 
and  f  requency  response  f  unctions  a  re  generated.  M  any  m  ethods  a  re  a  pplicable  for  m  odal 
identification.  However,  for  Mir,  the  input,  ships  docking,  were  not  measured  making  common 
practices  void.  T he  method  investigated  to  solve  the  challenges  involved  with  the  Mir  was  the 
inverse  system  identification  technique.  T  his  approach  initially  presented  by  Kammer  [1999] 
identifies  the  vibrational  characteristics  of  an  inverse  representation.  An  inverse  representation 
differs  from  a  normal  system  representation  in  that  the  input  and  output  are  switched.  Meaning, 
the  changes  in  the  natural  frequencies  and  modes  directly  relate  to  the  system  parameters  of  the 
structure,  for  a  more  detailed  look  refer  to  Kammer  et  al.  [1999].  A  2646  DOF  FE  model  of  Mir 
was  utilized  to  match  the  analytical  output  to  the  measured.  From  this  the  input  was  derived  and 
more  thorough  tests  on  the  structural  health  of  Mir  were  conducted. 

The  1  argest  i  nterplanetary  s  pacecraft  ev  er  d  eveloped,  C  assini,  1  aunched  o  n  a  T  itan-IV, 
Centaur  1  aunch  v  chicle  r  equired  a  t  est  ve  rified  f  inite  e  lement  m  odel  t  o  be  c  ompleted  and 
approved  by  NASA.  T  he  validation  of  the  model  was  critical  to  final  verification  of  coupled 
loads  analysis  and  margin  of  safety  assessments,  Coleman  et  al.  [1996].  Modal  data  was  taken 
using  a  s tatic  t  est  and  c  omponent  s  ine  s  weep  t  ests  t o  be  us  ed  t o  v  erify  t he  F  E  m  odel.  Before 
testing  could  occur,  an  FE  model  with  a  correct  mass  and  configuration  was  used  to  determine 
shaker  location  and  assess  instrumentation  location.  The  FE  model  Guyan-reduced  mass  matrix 
from  the  FE  model  was  used  for  orthogonality  and  effective  mass  calculations.  The  analytical 
approach  used  to  correct  the  FE  model  was  the  eigenvalue  problem,  this  method  allows  for  the 
variation  of  m  any  m  odal  pa  rameters  toe  nsure  a  ve  rified  m  odel.  T  he  obj  ective  f  unction 
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represented  b  elow  was  different  b  eeause  i  t  w  as  b  ased  p  rimarily  o  n  f  requencies;  t  hey  were 
stressed  to  be  the  most  important. 

«=  Z  '*',■.(//•//>"+  Z  <2) 

pairs  _j,k  pairs  _j,k 

where  Wfj.  and  are  diagonal  weighting  matrices  for  frequeney  and  eross  orthogonality, 

with  nonzero  values  at  matching  test  mode  pairs.  Also,  ff  and  //are  the  analytical  and  test 

frequencies,  Xj^  is  the  cross  orthogonality  matrix  between  analytical  modes  j  and  test  modes  k. 

Xj^  is  th  e  d  esired  c  ross  o  rthogonality  ma  trix.  O  ptimal  mo  del  p  arameters  were  achieved  by 

driving  this  error  function  e  down  to  zero  by  updating  model  idealizations,  joint  flexibilities  and 
stiffness. 

The  tuning  approach  that  m ost  cl osely  resembles  t he  m ethod u sed  to  tune  PnPSAT  I  was 
developed  by  Doupe  et  al.  [2009].  They  created  FE  models  of  the  panels  of  Falcon  SAT  5  (FS5), 
experimentally  collected  natural  frequencies  and  modes  using  a  scanning  laser  vibrometer  of  the 
impulse  response  of  each  panel.  Every  panel  FE  model  was  tuned  and  then  a  full  FS5  model  was 
created  b  y  connecting  t  he  p  anels  together  b  y  spring/damper  e  lements.  Next,  th  e  imp  ulse 
response  was  collected  for  the  fully  constructed  FS5  and  was  used  to  tune  the  completed  FS5  FE 
model.  The  ad  vantage  of  t  his  p  rocess  was  t  he  act  ual  s  atellite  an  d  al  1  i  ts  co  mponents  w  ere 
available  t  o  b  e  t  aken  a  part  an  d  t  ested.  T  herefore,  all  masses  could  be  m  easured  making 
modeling  the  mass  more  certain.  Another  advantage  of  possessing  the  entire  satellite  is  that  it  is 
possible  to  gather  not  only  the  natural  frequencies  but  also  the  lower  frequency  modes  of  the 
constructed  satellite.  This  allows  the  objective  function  to  minimize  the  difference  between  the 
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modes  w  hile  ma  intaining  fixed  natural  f  requency  co  nstraints  1  eading  t  o  an  accu  rately  t  uned 
complete  FE  model  for  0-300  Hz  for  five  modes. 

2.2  Finite  Element  Analysis 

FE  analysis  is  a  numerical  process  involving  solution  of  field  problems  Cook  et  al.  [2002], 
The  field  values  are  any  dependent  variables  described  by  differential  or  integral  equations.  FE 
codes  reformulate  d  ifferential  o  r  in  tegral  e  quations  in  to  algebraic  o  nes  using  i  nterpolation 
functions  and  solve  them  at  discrete  points  called  nodes.  Each  node  is  connected  to  other  nodes 
through  elements.  E  lements  c  an  b  e  composed  of  various  geometric  s  hapes  from  a  one 
dimensional  to  a  three-dimensional  element.  An  arrangement  of  elements  creates  a  mesh  which 
represents  the  structure.  In  creating  an  FE  model,  many  common  assumptions  are  made  to  avoid 
over  c  omplication  of  the  F  E  f  unctions  mesh.  F  irst,  c  omplex  geometric  shapes  like  isogrid 
panels,  hi  nges  a  nd  c  omponents  are  u  sually  s  implified  to  simpler  shapes  such  ass  quares  a  nd 
rectangles.  In  the  case  of  the  PnPSAT  I  model,  most  of  the  connectors  such  as  screws,  nuts  and 
bolts  ar  e  not  di  rectly  m  odeled  t  o  reduce  m  odel  c  omplexity.  Second,  m  aterials  ar  e  as  sumed 
homogenous,  m  eaning  t  hey  possess  t  he  s  ame  material  pr  operties  t  hroughout.  PnPSAT  I  is 
predominantly  an  all  aluminum  structure.  Third,  materials  are  assumed  isotropic  and  constant 
meaning  they  have  the  s  ame  properties  in  every  direction  and  their  properties  do  not  change. 
Small  displacements  and  rotations  are  assumed  along  with  fixed  loads  and  boundary  conditions. 
All  linear  assumptions  apply  with  respect  to  material  properties,  geometry,  and  loads. 

The  F E  an alysis  p rocess  generally  involves  three  steps:  pre-processing,  numerical  analysis, 
and  post-processing.  Pre-processing  involves  creating  an  FE  model  of  the  test  article  to  include 
choosing  e  lement  t  ypes,  m  aterial  pr  operties,  1  oads,  bounda  ry  c  onditions,  a  nd  m  esh  de  nsity. 
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Next,  an  eigenanalysis  is  completed  which  acquires  the  modes  and  natural  frequencies.  Finally, 
post-processing  is  performed  to  analyze  the  computed  analytical  results. 

2.3  Modal  Analysis 

Modal  an  alysis  i  s  used  to  measure  the  d  ynamic  r  espouse  of  P  nPSAT  1 1  o  va  rious  1  oads. 
Modal  an  alysis  i  s  t  he  p  rocess  o  f  m  easuring  t  he  t  hree  d  ynamic  ch  aracteristics  o  f  s  tructures: 
modes,  damping  ratios,  and  natural  frequencies.  A  mode  is  the  magnitude-independent  vector 
representation  of  the  shape  of  a  structure  when  the  structure  is  resonating  at  a  natural  frequency. 
Damping  involves  any  mechanism  that  removes  energy  from  a  structure.  F or  example,  a  shock 
absorber  is  a  classic  example  of  a  viscous  damper.  Natural  frequencies,  also  known  as  resonant 
frequencies,  a  re  t  he  f  requencies  i  n  w  hich  t  he  s  tructure  h  as  a  t  endency  too  scillate  at  1  arger 
amplitudes. 

The  discrete  FE  equations  of  motion  (EOM)  for  complex  structures  can  be  written  in  the  time 
domain  as: 


M  x{t)  +  {\  +  iY)Kx{t)  =  F{t)  (3) 

where  MeR™  represents  the  mass  matrix,  represents  the  stiffness  matrix  and  y  is 

the  structural  damping  coefficient.  Both  M  and  K  matrices  are  symmetric.  The  reduced  PnPSAT 
I  FE  model  developed  in  this  research  has  over  40,000  elements  and  70,000  nodes.  E  ach  node 
possesses  either  three  or  six  DOFs  making  n  »  200, 000 .  F{t)  represents  the  force  matrix  and  is 

F  e  R“'"  which  can  have  tens  to  hundreds  of  load  cases.  Because  PnPSAT  I  is  predominately 
aluminum,  onl  y  on  e  s  tructural  da  mping  coefficient  /  is  u  sed.  U  sually,  a  viscous  da  mping 
matrix  C  is  present  in  Eq  (1);  however,  there  are  no  viscous  dampers  in  PnPSAT  I,  damping  is 
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purely  a  f unction  of  t  he  c  omplex  m  odulus,  iy .  The  x{t)  e  and  x{t)  e  represent  t  he 
acceleration  and  displacement  vectors,  respectively  Dowpe  [2009], 

For  the  undamped  free  vibration  case  Eq.  (3)  can  be  written  as 

M  x{t)  +  Kx{t)  =  0  (4) 

First,  a  synchronous  solution  is  assumed  of  the  form, 

x{t)  =  (5) 

where  ®  is  the  radian  frequency  and  ^  is  a  co  nstant  real-valued  vector  of  dimension  n.  A  fter 
substitution,  a  real-valued  eigenvalue  problem  (EVP)  is  formulated 

K'^  =  (6) 

where  the  eigenvalue  A  =  o/ .  Eq.  (6)  can  be  solved  for  the  n  eigenvectors  (j)  for  n  values  of  A . 
However,  t  ypically  onl  y  a  few  oft  he  t  otal  num  ber  of  e  igenvectors  a  nd  c  orresponding 
eigenvalues  less  than  a  predetermined  cutoff  are  actually  computed.  Combining  the  m  computed 
eigenvectors  (j)  into  a  matrix  yields 

U  =  ...,  (7) 

Both  the  mass  and  stiffness  matrices  are  projected  onto  the  eigenspace  by 

K'  =  U^KU  (8) 

M'  =  U^MU  (9) 

The  now  diagonalized  mass  and  stiffness  matrix  allows  for  easier  computation  of  the  frequency 
response  of  the  structure  from  an  input  force  F{co)  from 
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X{ia))  =  U{-a/M '+  (1  +  iY)K  ')"*  U^F{a)). 


(10) 


The  magnitude  \X{ia>)\  and  phase  angle  of  ZX(io))  is  eomputed  from 

I  X(io})  \=  ^[RQ(X(io}))f  +  lm(X(io})f  (11) 


p(co)  =  tan 


-lm(X  (ico)) 
RQiX(ico)) 


(12) 


where  Re  extraets  the  real  values  and  Im  extraets  the  imaginary  values. 

Plots  of  magnitude  |  X(ico)  \  and  phase  ZX(ia>)  versus  co  are  useful  tools  in  modal  analysis. 
The  peaks  located  in  the  magnitude  plot  indicate  where  the  natural  frequencies  of  the  structure. 
Also,  the  modes  can  be  extracted  from  the  magnitude  and  phase  information.  In  an  optimization 
process,  t  he  e  igenvalues  a  nd  eigenvectors  of  t  he  F  E  m  odel  are  t  uned  i  n  or  der  tom  atch  t  he 
experimentally  measured  natural  frequencies  and  modes,  respectively  Dow/je  [2009]. 


2.4  Scanning  Laser  Vibrometer  Data  Collection  Fundamentals 

A  few  different  methods  can  be  used  to  acquire  frequency  response  function  (FRF)  from  a 
structure.  T  he  m  ost  c  ommon  m  ethod  i  nvolves  the  us  e  of  a  ccelerometers.  A  ccelerometers,  a  s 
their  name  suggests,  measures  the  acceleration  of  the  subject  being  tested.  A  ccelerometers  are 
necessary  in  providing  the  data  to  create  FRFs  at  each  location  of  which  have  to  be  combined  to 
create  the  overall  mode  shapes  of  the  test  subject. 

Typically,  dozens  or  more  accelerometers  are  placed  over  a  large  test  article.  The  addition  of 
accelerometers  can  change  the  response  of  the  system  depending  on  the  mass  and  stiffness  of  the 
test  subject.  H  owever,  in  order  to  acquire  an  accurate  representation  of  the  mode  shapes  using 
accelerometers,  hundreds  of  accelerometers  may  be  required  depending  on  the  complexity  of  the 
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structure  and  accuracy  required.  This  is  why  modem  modal  analysis  is  turning  to  the  use  of  laser 
Doppler  vibrometer  (LDV). 

LDV  p  rovides  m  any  d  ifferent  b  enefits,  of  w  hich  m  easuring  test  ar  tides  without  contact 
provides  the  1  argest  b  enefit.  Also,  i  t  i  s  po  ssible  t  o  easily  measure  a  t  tens  tot  housands  of 
locations  enabling  a  more  accurate  extraction  of  the  modes,  especially  higher  frequency  modes. 

The  Polytech  scanning  laser  vibrometer  (PSV),  model  PSV-400-3d-M,  is  shown  in  Figure  6, 
and  was  used  to  collect  data  for  modal  analysis.  The  three  laser  heads  are  focused  on  the  same 
point  a  nd  ut  ilize  a  D  oppler  e  ffect  m  ethod  f  or  measuring  v  elocity.  T  he  D  oppler  e  ffect  i  s  t  he 
change  in  frequency  of  a  wave  for  an  observer  moving  relative  to  the  source  of  the  wave.  The 
most  c  ommon  e  xample  of  this  phenomenon  i  s  when  a  ve  hide  s  ounding  i  ts  s  iren  approaches. 
The  s  iren  w  ill  be  hi  gher  pi  tched  w  hen  t  raveling  t  owards  t  he  obs  erver  and  1  ower  pi  tched  a  s  i  t 
travels  away.  U  tilizing  the  D oppler  e quations,  the  vi brometer  e mits  laser  energy  at  a  known 
frequency  and  m  easures  t  he  f  requency  of  t  he  r  eflected  1  ight  t  o  de  termine  th  e  v  elocity  o  f  th  e 
particular  scan  point  on  the  test  article  Doupe  [2009].  The  Doppler  equation,  adopted  from  Rees 
[2001]  is 


^  vcos6> 


c 


(13) 


where  c  is  the  speed  of  light,  v  is  the  velocity  of  the  scan  point,  0  is  the  angle  of  the  detector  in 
relation  to  the  source,  is  the  transmitted  frequency,  X^  is  the  detected  frequency. 
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In  the  testing  of  PnPSAT  I  panels  it  was  necessary  only  to  use  one  laser  head  to  measure  the 
out-of-plane  displacements.  M  ensuring  the  in-plane  displacements  proved  unnecessary  because 
the  in-plane  modes  have  natural  frequencies  that  far  exceed  the  range  of  tuning  conducted  in  this 
research.  After  complete  alignment  of  the  laser  head,  a  grid  of  known  coordinates  from  the  FE 
model  are  loaded  into  the  PSV  software  predetermining  the  measurement  locations.  Detailed 
information  on  the  Polytec  3D  LDV  and  PSV  software  can  be  found  in  Polytec  [2007a]  Polytec 
[2007c]. 

The  acquired  time  domain  data  from  the  scanning  laser  vibrometer  must  first  be  converted  to 
the  f  requency  dom  ain.  T  he  PSV  software  co  nverts  t  he  t  ime  r  esponse  d  ata  to  the  c  omplex 
frequency  domain  through  a  fast  Fourier  transform  (FFT).  The  continuous  FFT  is  represented  by 


F{a))  =  [  “  fit)e~‘^‘dt 

J  -00 


(14) 
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where  F{co)  is  t  he  f  requency  dom  ain  s  ignal,  / {t)  is  t  he  t  ime  dom  ain  s  ignal,  a  nd  ft)  is  th  e 
frequency.  H  owever,  s  ince  d  iscrete  s  ampling  t  imes  are  used  t  o  acq  uire  t  he  d  ata  t  he  d  iscrete 
Fourier  t  ransform  (  DFT)  m  ust  be  us  ed.  T  he  D  FT  c  onverts  di  screte  t  ime  dom  ain  da  ta  i  nto 
discrete  frequency  domain  data.  By 

F„="^f{t,)e-^"\n  =  Q...,N-\  (15) 

i=0 

where  N  is  the  number  of  samples.  In  order  to  prevent  aliasing,  which  occurs  when  two  signals 
of  d  ifferent  fre  quencies  are  s  ampled  at  s  uch  a  r  ate  as  t  hey  are  recorded  a  s  t  he  s  ame  s  ignal,  an 
anti-aliasing  f  ilter  i  s  us  ed.  T  he  P  SV  au  tomatically  s  ets  t  he  co  rrect  an  ti-alias  f  ilter  o  nee  th  e 
desired  frequency  range  is  established  by  the  user. 

Leakage  is  another  cause  for  error  and  is  an  effect  of  signals  which  are  not  periodic  in  the 
sampled  data.  It  is  caused  by  the  finite  of  the  data  record,  that  introduce  discontinuities  across 
stop/start  regions.  In  order  to  remedy  leakage,  a  window  is  applied  to  the  time  response  data. 
Various  windows  have  been  created  for  this  purpose  to  include  rectangular,  Hanning,  Kaiser- 
Bessel,  flat  top,  and  exponential.  T  he  windows  of  choice  in  impact  testing  are  the  force  and 
exponential  w  indows.  T  he  f  orce  w  indow  c  atches  t  he  qui  ck  i  mpulse  of  t  he  s  triking  ha  mmer. 
While  the  exponential  window  captures  the  output  which  is  a  large  magnitude  sinusoid  initially 
because  i  t  w  as  s  truck  b  y  t  he  h  ammer  t  hen  i  t  ‘  exponentially’ t  apers  o  ff  t  owards  t  he  end  of  t  he 
sample. 

The  complex  frequency  data  is  converted  into  an  FRF.  W ith  the  frequency  domain  data  the 
vibrometer  software  calculates  the  FRF  using 

FRF  =  °‘“P‘“FFT  (IS) 

inputFFT 
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Unfortunately,  this  computation  is  more  difficult  than  it  appears.  First  the  correlation  functions 
must  be  determined  using  the  general  form 


\  f(T^)s{t  +  '^)dT  (17) 


from  this  the  power  spectral  densities  can  be  found  by  doing  the  FFT  of  the  correlation  functions. 
The  four  power  spectral  densities  with  respect  to  the  input  R{co)  and  the  output  C{co)  are 

Scc{o))  =  C{(o)C{co) 

S„r{(o)  =  R{(o)R((o) 

—  Uo) 

S^i^{co)  =  C{co)R{co) 

S^^{co)  =  R{co)C{co) 

were  R{co)  and  C{co)  are  the  FFTs  of  the  input  and  output  respectively.  The  over  bar  represents 
the  complex  conjugate.  Given  these  power  spectral  densities,  it  is  possible  to  get  two  FRFs,  HI 
and  H2  given  by 


Hl{o))  = 


^CR  (^) 
^RR 


and  H2ico)  = 

^RC  (®) 


(19) 


Both  FRFs  represent  some  form  of  output  divided  by  input,  but  the  HI  FRF  is  primarily  affected 
by  noise  on  the  input.  The  H2  estimator  is  primarily  affected  by  noise  on  the  output,  hr  general 
output  noise  is  greater  than  input  noise;  hence,  the  H\  FRF  is  better  for  viewing.  However,  both 
FRFs  should  be  almost  identical  when  plotted  so  an  averaging  technique  is  used  to  minimize  the 
noise.  Each  scan  point  is  measured  ten  times;  consequently  each  scan  has  four  power  spectral 
densities.  Each  power  spectral  density  is  averaged  with  the  others  of  its  same  variable  name  then 
the  averaged  FRFs  HI  and  H2  are  computed. 
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In  order  to  ensure  that  the  data  is  aceurate,  a  test  parameter,  eoherence,  is  used.  Coherence  is 
the  ratio  of  H1/H2  and  is  a  measure  of  consistency  in  the  measured  data.  It  relates  how  much  of 
the  output  is  directly  caused  by  the  input.  Coherence  ranges  from  zero  to  one,  one  being  the 
input  and  output  are  perfectly  related  and  zero  being  the  input  and  output  are  totally  unrelated. 

After  the  scanning  LDV  has  been  setup  correctly,  it  can  be  put  in  a  scan  mode  where  it  takes 
ten  measurements  at  each  scan  point  until  completion.  The  PSV  software  will  classify  the  data 
based  on  valid,  optimal,  low  signal,  over  range  or  invalidated.  The  over  range  data  appears  when 
the  magnitude  of  the  velocity  is  greater  than  the  selected  laser  sensitivity.  It  is  optimal  to  run  the 
testing  at  night  in  order  to  minimize  the  effects  of  ambient  noise.  The  scanning  laser  vibrometer 
and  test  article  are  very  susceptible  to  noise  to  include  anything  from  people  walking  near  it  or 
talking  in  the  s  ame  room.  O  nee  completed  the  data  i  s  processed  and  u  sed  for  tuning  the  FE 
model. 

2.5  Finite  Element  Model  Tuning 

Costly  and  time  consuming  AIT  testing  can  possibly  be  reduced  or  even  avoided  by  the  use  of 
accurately  tuned  FE  models.  Tuning  essentially  means  a  model  is  corrected  to  accurately  predict 
experimental  data  over  a  range  o  f  excitation  frequencies.  A  tuned  model  should  theoretically 
allow  for  the  prediction  of  natural  frequencies  and  modes  of  various  configurations.  Creating  an 
FE  model  can  vary  significantly  in  difficulty.  Regardless  of  complexity,  creating  an  accurate 
model  on  the  first  attempt  is  hardly  if  ever  realized.  The  process  of  tuning  adjusts  parameters  in 
the  FE  m  odel  to  m  ore  a  ccurately  pr  edict  r  esponses  t  o  va  rious  1  oads.  Tuning  c  an  be 
accomplished  s  everal  d  ifferent  w  ays  t  o  i  nclude  c  hanging  parameters  which  af  feet  t  he  m  ass, 
stiffness,  or  damping  matrices.  More  commonly,  the  modulus  of  elasticity  or  material  density;  is 
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altered.  In  this  research,  only  the  modulus  of  elasticity’s  of  the  adapter  ring  and  spring  stiffness 
coefficients  between  the  baseplate  and  -Z  panel  will  be  altered.  A  brief  overview  of  FE  model 
tuning  is  presented,  next  while  the  details  can  be  found  in  Chapter  III. 

In  o rder  t o  ac curately  t une  the  FE  m odel,  experimental  m odal  d ata  o f  t he  test  a rticle  f irst 
must  be  co  llected.  To  a  cquire  t  he  m  odes  o  f  t  he  s  tructure  e  xperimentally,  the  boundary 
conditions  must  be  known.  Boundary  conditions  typically  support  or  constrain  the  structure.  An 
example  of  a  constrained  structure  is  a  satellite  fixed  on  a  shaker  table  to  simulate  the  attachment 
to  a  rocket.  In  certain  circumstances,  it  is  better  to  allow  the  structure  to  float  in  order  to  keep 
the  dynamic  properties  independent  of  any  constraints  or  supports. 

There  are  several  devices  used  to  excite  a  structure  to  include  shaker  tables,  electromagnetic 
shakers,  impact  hammers,  just  to  name  a  f  ew.  Modes  whose  natural  frequencies  are  near  the 
excitation  frequency  are  excited.  O  nee  the  structure  is  excited  it  is  important  to  measure  the 
input  excitation  and  out  put  response.  As  m  entioned  e  arlier,  t  he  out  put  can  b  e  measured  b  y 
accelerometers  or  laser  vibrometers.  Now  an  input-output  relationship  can  be  formulated  in  the 
frequency  domain  via  FRFs  where  the  peaks  of  the  FRFs  are  the  natural  frequencies,  as  shown  in 
Figure  7. 


Vib  Top  &  Reft  HI  Accele. . . 
Frequency 

100.0  Hz 

■  Magnitude 

10.293  dB 


Figure  7:  Frequency  Response  Function  of-Z  panel 


23 


Damping  can  al  so  b  e  c  alculated  from  t  he  F  RFs  using  t  he  ha  If-power  bandwidth  method. 
Taking  me  asurements  a  t  mu  Itiple  locations  a  Hows  f  or  th  e  c  reation  o  f  operational  de  flection 
shape  (ODS).  S  oftware  utilizing  a  curve-fitting  approach  allows  for  the  extraction  of  modes, 
damping  factors,  and  natural  frequencies  from  the  FRFs  and  ODSs. 

Upon  extraction  of  the  experimental  data,  tuning  of  the  FE  model  can  begin.  The  concept  is  to 
modify  the  eigenvectors  to  match  the  measured  modes  and  the  eigenvalues  to  match  the  square 
of  t  he  n  atural  frequencies.  In  t  he  t  uning  pr  ocess,  parameters  af  fecting  m  ass  a  re  not  altered 
because  the  mass  is  an  easily  measurable  quantity  with  a  high  degree  of  accuracy.  Both  of  these 
parameters  affect  the  stiffness  matrix  which  in  turn  changes  the  eigenvectors  and  eigenvalues. 

Due  to  the  complexity  and  size  of  FE  models  it  is  nearly  impossible  to  tune  the  FE  model  by 
hand.  It  i  s  c  ommon  t  o  minimize  an  obj  ective  f  unction  which  i  s  a  m  easure  of  t  he  di  fference 
between  m  easured  an  d  analytical  r  esponses.  A  c  ommon  obj  ective  function  i  s  as  previously 
presented  and  repeated  here  the  difference  between  the  measured  and  analytical  eigenvalues  and 
eigenvectors  given  as 


J  =  'La, 


1- 


P  r 

/=i  j=\ 


1- 


d).. 

n/ 


(1) 


The  weighting  factors  a,  and  bj  can  be  adjusted  to  create  a  variety  of  objective  functions  for 
methods  of  comparison  that  place  more  importance  on  a  particular  mode  or  natural  frequency, 
and  are  the  i"'  analytical  and  measured  natural  frequencies,  respectively,  and  and 

are  t  he  i‘’'  analytical  an  d  m  easured  e  igenvectors  or  m  odes,  r  espectively.  T  he  an  alytical 

A 

eigenvectors  are  normalized  with  .  The  scalar  summation  limits  p  and  r  are  the  number  of 
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desired  modes  and  desired  number  of  measurement  locations,  respectively.  The  absolute  value 
operation  on  t  he  1  ast  t  erm  i  s  r  equired  be  cause  t  he  overall  s  ign  o  f  a  p  articular  ei  genvector  i  s 
unknown.  The  objective  function  J  is  minimized  changing  the  parameters  chosen  while  obeying 
the  constraints  in  the  FE  model.  By  minimizing  the  objective  function,  the  differences  between 
analytical  and  m  easured  natural  f  requencies  a  nd  m  odes  a  re  r  educed.  This  m  inimization 
approach  can  match  natural  frequencies  quite  easily,  because  they  are  simply  scalar  values,  but  it 
is  much  more  difficult  to  minimize  the  differences  between  modes  which  are  vector  quantities. 
In  t  he  ne  xt  c  hapter  t  he  de  tails  i  nvolving  t  esting,  m  odeling,  a  nd  t  uning  w  ill  be  pr  esented. 
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III.  Method 


This  chapter  first  describes  the  methods  used  to  generate  an  untuned  FE  model  for  the  -Y  and 
-Z  p  anels.  N  ext,  t  his  ch  apter  d  iscusses  h  ow  ex  perimental  d  ata  is  collected  tot  une  p  anel  FE 
models.  Finally,  the  PnPSAT  I  FE  model  is  tuned  and  tested. 

3.1  Method  Overview 

In  order  to  create  an  FE  model  of  PnPSAT  I,  only  FE  models  of  the  -Y  and  -Z  panels  are 
created.  Modal  testing  is  conducted  to  collect  modes  and  natural  frequencies  for  both  the  -Y  and 
-Z  panels.  The  description  of  this  testing  will  include  the  equipment  used,  test  setup  and  theory. 
Also,  the  tuning  process  was  applied  to  the  -Y  and  -Z  panels  is  described  in  detail.  The  tuned 
panel  m  odels  are  c  opied  an  d  at  tached  t  o  cr  eate  a  structure  Pn  PSAT  I  FE  m  odel.  S  ix  DOF 
spring/damper  elements  keep  the  base  plate  and  the  adapter  ring  attached.  T  he  -Z  panel,  base 
plate,  and  component  p  oint  m asses  are  c  onnected  vi  a  rigid  1  inks  p roviding  a  fully  de  veloped 
PnPSAT  I  FE  model.  This  FE  model  will  be  tuned  once  using  this  configuration  to  match  KAFB 
experimentally  collected  data  from  the  ORS  office. 

3.2  Generate -Y  and  -Z  Panels  FE  Models 

Modeling  i  s  s  een  as  a  way  tor  educe  c  omplex  a  nd  expensive  t  esting.  P  revious  t  o  t  his 
research  effort  a  high  order  FE  model  was  created  by  Spaceworks.  High-order  models  are  easily 
created  from  a  detailed  CAD  file  using  automatic  mesh  generation  f eature  f ound  i n  m ost  F E 
packages  today.  These  auto-generated  meshes  generally  possess  many  more  nodes  and  elements 
than  those  that  are  carefully  controlled  by  the  analyst.  Although  generating  a  high-order  model 
is  initially  less  time  consuming,  the  fact  that  the  Spaceworks  FE  model  has  more  than  a  million 
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DOF  results  i  n  hi  gh  de  mands  on  computing  pow  er  t  o  s  olve  t  he  e  igenvalue  pr  oblem.  F  or 
example,  approximately  18  m  inutes  on  4  dua  1  core  AMD  8200  2.8G  Hz  pr  ocessors  w  ith  64 
gigabytes  of  RAM  and  requires  9.634  gigabytes  of  scratch  disc  space.  On  the  same  machine  an 
eigenanalysis  of  the  PnPSAT  I  can  be  completed  in  3  minutes.  The  factor  of  six  difference  might 
not  s  eem  s  ignificant  but  dur  ing  t  he  opt  imization  pr  ocess,  this  e  igenvalues  a  nd  nor  mal  m  odes 
analysis  is  calculated  hundreds  to  thousands,  so  of  times  this  factor  makes  a  significant  impact. 
Also,  during  the  optimization  process,  files  are  read  and  saved  to  disc  because  the  matrices  are 
too  large  to  reside  in  RAM  which  dramatically  slows  down  the  tuning  process  because  disc  read 
write  speed  is  dependent  solely  on  the  hard  disc.  Therefore,  a  reduced  order  FE  model  needs  to 
be  created  for  the  tuning  process  to  be  computed  in  a  reasonable  amount  of  time. 

The  first  step  is  to  measure  the  mass  of  the  -Y  and  -Z  panels.  When  the  FE  model’s  mass  is 
in  close  agreement  with  the  measured  mass  then  only  the  stiffness  matrix  has  to  be  tuned.  The 
+X,-X,-Y,  +Y  panels  are  constructed  in  similar  ways  possessing  the  same  properties  to  include 
structure,  electronics,  mass.  The  -Z,  +Z  are  constructed  similarly.  Hence,  tuning  the  -Y  and  -Z 
panels  is  sufficient  because  once  the  two  panels  are  tuned,  they  can  be  recopied  and  combined  to 
create  the  outer  structure  of  PnPSAT  I.  A  ddition  of  the  components  via  rigid  links  allows  for 
accurate  t  uning  for  t  he  completed  P  nPSAT  I .  An  ove  rview  of  t  he  pr  ocess  i  s  i  n  Figure  5,  i  n 
Chapter  I. 

Each  PnPSAT  I  p  anel  is  b  ased  o  n  a  cl  amshell  design.  The  e  lectronics  a  re  hous  ed  i  n  t  he 
middle  and  the  two  outer  panels  are  fastened  together. 
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Unique 

Ribbed 

Straeture 


Figure  8:  Picture  of  one  half  of  PnPSAT  I  -Y  panel  with  electrical  components;  arrows  point  out 
the  unique  ribbed  structure  which  joins  the  two  halves  of  the  panel  together. 


Notice  the  ribbed  interior  and  exterior  of  the  panel.  This  design  is  unique  in  that  it  protects 
the  electrical  components  and  results  in  a  very  stiff  but  lightweight  panel.  C  reation  of  the  FE 
model  of  PnPSAT  was  initially  generated  from  a  detailed  CAD  file  provided  by  Spaceworks. 
From  the  CAD  model  of  -Y  and  -Z  panels,  reduced  FE  models  were  created.  F  EMAP,  is  a 
useful  t  ool  i  n  creating  a  nd  analyzing  f  inite  element  m  odels  t  hrough  a  us  er  friendly  GUI 
(graphical  user  interface).  The  top  shell  of  the  -Z  panel  consists  of  geometry  points,  surfaces 
and  solids  shown  in  Figure  9. 


28 


Figure  9:  CAD  geometry  model  of  the  -Z  panel. 


The  next  task  consisted  of  creating  a  geometrically  accurate  FE  model  from  the  CAD  files. 
Nodes  are  generated  individually  to  ensure  a  continuous  and  accurate  mesh.  S  tailing  with  the 
comers  a  nd  s  ignificant  pa  its  of  t  he  ge  ometry  t  he  node  s  a  re  pi  aced  i  n  a  w  ay  as  t  o  e  nsure 
continuity  and  dynamics  accuracy.  Each  node  in  Figure  10  was  deliberately  placed  in  order  to 
ensure  e  lement’s  node  s  w ould  be  c  oincident  toe  ach  ot her.  This  en  sures  el  ements  ar  e  r  igidly 
attached  and  are  not  separating  from  one  another. 
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Figure  10:  Nodal  representation  of  the  -Z  panel. 
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Because  the  majority  of  PnPSAT  I  is  mo  deled  as  6061-T651  a  luminum,  each  dementis 
modeled  w  ith  th  e  s  ame  ma  terial  properties;  h  owever,  over  60  material  property  v  alues  are 
included  in  order  to  allow  the  tuning  algorithm  to  adjust  the  stiffness  properties  of  individual 
sections  of  t  he  pa  nel  m  odds.  Therefore,  the  Young’s  m  odulii  can  be  a  djusted  i  n  t  he 
optimization  p  rocess  to  correlate  th  e  F  E  mo  del  and  th  e  e  xperimentally  collected  modes  an  d 
natural  frequencies. 


Figure  11:  -Z  Panel  FE  Model  where,  each  color  represents  a  different  material  card  (95  total); 
these  material  card’s  properties  are  altered  in  the  tuning  process. 


Approximately  89%  of  the  elements  are  solid  8  node ,  24  D  OF  hex  elements.  Each  node  of 
the  h  ex  s  olid  e  lement  h  as  3  tr  anslational  D  OFs.  It  is  e  asy  to  generate  a  F  E  m  esh  o  f  s  olid 
elements  from  already  created  nodes.  This  generation  is  simply  connecting  the  dots  on  a  large 
scale.  Approximately  10%  of  the  elements  are  3  node  18  DOF  plate  elements,  which  are  used  on 
the  exterior  structure  of  the  panel.  Ensuring  the  thickness  of  these  plate  elements  agree  closely 
with  the  actual  satellite  panel  thickness  was  critical  in  matching  the  mass  and  more  importantly, 
that  the  m  ode  1  ocations  are  i  n  cl  ose  a  greement  w  ith  t  he  ex  perimentally  measured  d  ata.  T  he 
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remaining  elements  w  ere  m  ade  up  of  r  igid  1  ink  e  lements  w  hieh  s  eeured  t  he  hi  nges  t  o  t  he 
structure  of  the  satellite. 

A  problem  arose  in  the  computation  of  the  FE  models  mass  matrix  because  the  panel  was 
not  constrained.  The  M SC  N astran  algorithm  that  creates  the  mass  matrix  produces  a  singular 
mass  matrix  when  it  s  houldn’t.  T  his  created  havoc  on  calculation  of  the  stiffness  matrix  and 
tuning  process  because  it  ruins  the  cross-orthogonality  check  (as  is  explained  later).  To  remedy 
this,  spring  elements  were  inserted  into  the  comers  of  the  stmcture  that  have  very  low  spring 
stiffness  c  oefficients.  These  s  prings  c  hanged  t he  first  s  ix  r igid  bod  y  m  odes  t  o  flexible  m  odes 
with  very  small  natural  frequencies,  six  rigid  body  modes  between  6  and  10  Hz.  This  alteration 
slightly  changed  the  natural  frequencies  but  left  the  mode  shapes  in  the  correct  order. 


Figure  12:  Transparent  -Z  panel  showing  the  spring  dampers  (black)  in  the  comers  with  fixed 
node  constraint.  Elevates  rigid  body  mode  to  avoid  singular  mass  and  stiffness  matrices. 

Both  the  -Y  and  -Z  panels  were  constracted  similarly  and  the  next  step  is  to  tune  the  panel 
FE  models  but  first  the  experimental  data  must  be  collected. 
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3.S  Test  Setup 


The  PnPSAT  I  modal  testing  at  A  FIT  consists  of  using  a  scanning  LDV  to  collect  velocity 
measurements  over  a  gr  id  of  measurement  1  ocations  o  n  t  he  s  urface  o  f  t  he  s  atellite  p  anels. 
PnPSAT  I’s  panels  are  mechanically  coarse  structures  which  made  choosing  ideal  measurement 
points  difficult.  An  automated  impact  hammer  is  used  to  excite  as  many  vibrational  modes  as 
possible  in  t  he  s  pacecraft  pa  nels.  T  he  i  mpulsive  f  orce  i  mparted  t  o  t  he  P  nPSAT  I  pa  nels  i  s 
measured  an  d  used  as  t  he  r  eference  s  ignal  from  w  hich  FRFs  are  c  omputed  from  each 
measurement  point,  of  which  there  are  over  150  measurement  points  for  each  panel. 

The  test  harness  shown  in  Figure  13  is  designed  to  simulate  free  vibration  by  minimizing  the 
amount  of  strain  in  the  panel  boundary  conditions.  For  true  free-free  boundary  conditions,  the 
panel  ideally  should  be  allowed  to  float  freely  in  space.  This  is  accomplished  by  suspending  the 
panel  by  bungee  cords.  Two  aluminum  beams  with  pyramid  shaped  foam  blocks  on  top,  press 
on  the  bottom  of  the  satellite  providing  light  damping.  The  foam  pads  are  necessary  to  damp  out 
the  s  ystem  b  etween  s  trikes  f  rom  t  he  s  haker.  Elevating  t  he  f  oam  p  ads  p  rovides  i  ncreased 
damping.  However,  the  pads  are  only  adjusted  to  provide  the  minimal  amount  of  damping  to 
ensure  the  structure  has  stopped  vibrating  due  to  the  excitation  applied  by  the  hammer. 
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Figure  13:  Test  harness. 


Once  the  test  harness  was  setup,  each  panel  was  placed  on  the  bungee  cord  and  foam  pads.  A 
coating  of  a  light  scattering  material  is  typically  required  as  a  surface  treatment  when  collecting 
LDV  data;  however,  due  to  PnPSAT  being  flight  hardware  with  onboard  electronics,  applying  a 
coating  on  the  panels  was  not  an  option.  Next,  a  reference  coordinate  frame  for  the  scanning 
LDV  i  s  d  efmed.  From  the  FE  model,  approximately  150  node  s  a  re  selected  at  which 
measurements  are  taken.  The  selected  nodes  are  all  on  the  same  plane  and  are  located  at  least 
1cm  f  rom  t  he  e  dge  or  any  screw  h  oles  t  o  en  sure  t  hat  t  he  1  aser  h  ead  can  t  ake  u  nobstructed 
measurements.  The  next  step  is  a  careful  2D  alignment  of  the  laser.  The  3D  coordinates  of  each 
selected  node  are  imported  into  the  PSV  software.  Only  one  of  the  three  laser  heads  are  used 
because  the  v elocities  are  only  measured  out-of-plane  making  the  other  two  vibrometer  heads 
unnecessary.  In  or  der  to  t  each  t  he  1  aser  v  ibrometer  t  he  t  est  ar  tide  co  ordinate  s  ystem,  several 
points  of  known  coordinates  are  required.  To  achieve  the  highest  accuracy  of  these  locating 
points,  precision  cal  ipers  were  u  sed  tom  easure  t  he  surface  location  of  f  our  node  s  ne  ar  t  he 
comers  of  the  panels.  Reflective  tape  and  a  pen  point  are  used  to  mark  the  location  of  these  four 
nodes  as  indicated  by  the  circles  drawn  in  Figure  14. 
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Figure  14:  PnPSAT  -Y  panel  with  precise  nodes  and  arbitrary  points,  reflective  tape  located  on 
comers. 

An  el ectromagnetic  s haker  set  to  act  as  an  automatic  ping  hammer  (Figure  15)  provides 
impulse  like  excitiations.  A  Hewlett  Packard  33120A  15MHz/Arbitray  waveform  generator  is 
programmed  to  generate  a  narrow  pulse  given  a  specific  amplitude  and  frequency  width. 


Force  Cell 
Stinger 


Shaker 


Figure  15:  MB  Dynamics  Cal50  exciter  electroynamic  shaker,  setup  to  perform  as  an  automatic 
ping  hammer. 

These  impulsive  s  ingle  strike  impacts  are  s  et  t  o  occur  ev  cry  t  en  s  econds  to  allow  f  or  t  he 
system  to  damp  out  before  striking  again  and  theoretically  excite  all  frequencies.  However,  most 
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of  the  energy  is  shown  in  the  lower  frequeneies.  Impact  excitation  was  chosen  for  this  research 
primarily  because  the  excitation  device  is  not  attached  to  the  test  article  which  would  change  the 
response  o  f  t  he  s  tructure.  Impact  excitation  al  so  provides  m  uch  be  tter  c  oherence  t  han  ot  her 
methods .  As  pr  eviously  de  scribed,  coherence  i  s  t  he  m  easure  o  f  t  he  1  inear  co  rrelation  b  etween 
input  a  nd  out  put  or  ho  w  m  uch  t  he  r  espouse  o  f  a  s  tructure  de  pends  di  rectly  on  t  he  e  xcitation 
input.  A  display  of  the  coherence  during  the  test  of  the  -Z  panel  is  in  Figure  16. 


Figure  16:  Coherence  of  testing  on  the  -Z  panel. 

A  threaded  rod,  called  the  stinger,  connects  the  shakers  actuator  to  a  force  cell  which  has  a 
flat  plate  at  the  other  end  (Figure  15).  Transducers  inside  the  force  cell  vary  an  applied  voltage 
corresponding  to  t  he  i  mpact  f  orce  m  agnitude.  The  f  orce  cell  d  ata  i  s  s  ent  t  hrough  a  s  ignal 
conditioner,  which  amplifies  the  signal,  and  then  to  the  vibrometer  as  a  reference  input. 

The  shaker  is  placed  underneath  the  test  stand  oriented  so  that  it  strikes  the  panel  at  a  desired 
excitation  location  as  shown  in  Figure  17.  Multiple  strike  locations  are  necessary  to  ensure  that 
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all  the  modes  are  excited;  if  the  shaker  strikes  the  panel  at  a  modes’  nodal  line  that  mode  will  not 
be  seen  in  the  measured  deflection  shapes. 


(a)  (b) 

Figure  17:  (a)  Shaker  striking  points  on  bottom  side  of  panel  (tested  facing  other  direction)  (b) 
picture  of  the  fifth  mode,  dark  regions  represent  nodal  lines. 


An  important  parameter  to  monitor  on  the  input  excitation  signal  is  spectral  density.  This  is 
the  a  mount  of  pow  er  i  mparted  t  o  t  he  s  tructure  ate  ach  f  requency.  Friswell  and  Mottershead 
[1995]  provides  a  good  discussion  on  spectral  density  (aka  power  spectrum).  A  typical  spectral 
density  plot  is  shown  in  Figure  18.  A  general  rule  of  thumb  is  to  keep  the  spectral  density  within 
10  dB  to  20  dB  of  the  starting  value  throughout  the  frequency  of  interest.  If  the  spectral  density 
drops  of  f  t  oo  qui  ckly,  sufficient  e  nergy  w  ill  not  be  i  mparted  i  nto  t  he  s  tructure  at  hi  gher 
frequencies  r  esulting  i  n  p  oor  m  easurements  at  these  frequencies.  Spectral  d  ensity  c  an  b  e 
affected  by  hardness  of  the  impact  plate  connected  to  the  force  cell,  magnitude  and  duration  of 
the  impulse  from  the  shaker,  and  mass  ratio  of  test  article. 
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Figure  18:  Spectral  density  plot  of  an  impact  excitation  Doupe  [2009], 

Yet  another  important  characteristic  to  monitor  during  impact  excitation  is  the  time  domain 
signal  of  t  he  i  nput.  T  he  t  ime  dom  ain  s  ignal  s  hould  1  ook  a  s  c  lose  t  o  a  theoretical  i  mpulse  a  s 
possible.  However,  in  reality,  the  force  cell  will  output  small  changes  in  voltage  at  the  beginning 
and  e  nd  of  an  i  mpulse  strike  a  s  t  he  s  triker  be  gins  t  o  a  ccelerate  a  nd  again  a  s  i  t  r  etums  t  o  i  ts 
original  position  after  the  strike.  These  changes  appear  as  side  lobes  to  the  main  impulse  in  the 
time  domain  and  must  be  minimized  to  ensure  a  clean  strike.  The  anomalies  can  be  minimized 
by  adjusting  the  magnitude  and  frequency  settings  of  the  waveform  generator.  Large  side  lobes 
on  the  time  domain  input  have  an  adverse  affect  on  coherence  and  reduce  the  quality  of  the  data 
collected.  Adjustments  attempting  to  eliminate  side  lobes  can  also  lead  to  the  striker  hitting  the 
panel  multiple  times  during  one  impulse  because  the  impulse  width  is  too  large. 
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Figure  19:  Top  Left:  A  good  impact  excitation.  Top  Right:  A  squared-off  peak  on  pulse  indicates 
that  the  hammer  struck  the  test  article  too  slowly.  An  adjustment  of  the  pulse  width  is  necessary. 
Bottom  Left:  The  striker  hits  the  panel  multiple  times  indicating  both  frequency  and  magnitude 
needs  to  be  adjusted.  Bottom  Right:  S  ide  lobes  indicate  that  the  magnitude  is  too  high  Doupe. 
[2009] 
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The  p  anel’s  dynamic  r  esponse  is  c  ollected  over  a  r  ange  o  f  frequencies  from  0  -1000  Hz. 
Figure  20  plots  t  he  grid  of  m  easurement  1  ocations  that  w  ere  s  elected  in  the  F  E  m  odd  and 
exported  to  the  PSV  software  for  each  panel.  The  top  laser  head  of  the  vibrometer  focuses  on 
each  point  in  succession  and  records  the  velocity  as  a  function  of  time  as  the  panel  is  excited. 


Figure  20:  Typical  node  points  on  the  -Z  panel  from  P  SV.  Points  picked  from  FE  model  to 
avoid  holes  and  scratches.  Blue  means  unmeasured,  green  means  successful  measurement  and 
gray  means  over  ranging  typically  caused  due  to  people  talking  or  walking  around  the  room. 


The  excitation  force  for  each  impact  is  recorded  as  a  function  of  time  as  well.  Each  impact 
of  the  ping  hammer  triggers  the  vibrometer  to  collect  data  for  the  programmed  period  of  time  (10 
seconds)  and  with  the  programmed  sampling  rate.  The  sampling  rate  must  be  2x  greater  than  the 
signal  to  avoid  aliasing.  The  default  setting  in  the  PSV  software  for  the  sampling  rate  is  2.56 
times  the  maximum  frequency  of  interest  to  avoid  aliasing. 

The  test  setup  is  extremely  sensitive  in  that  it  cannot  be  touched  once  aligned  and  calibrated. 
The  ED  V  can  co  licet  v  ery  small  mo  tions,  f  or  example,  t  alking  a  nd  c  losing  door  s  are  eas  ily 
detected  see  Figure  20.  Due  to  sensitivity  to  noise,  tests  were  conducted  at  night  to  avoid  any 
erroneous  measurements.  Testing  at  n  ight  ensured  m  ore  accurate  r  esults  a  nd  a  voided  ove  r 
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ranging  of  the  lasers.  Each  point  was  measured  ten  times  and  averaged  in  the  frequency  domain 
to  min  imize  n  oise  e  ffects.  A  picture  of  th  e  c  omplete  te  st  s  etup  u  sing  th  e  -Y  p  and  shown  in 
Figure  21(a). 


Figure  21  (a)  Satellite  panel  in  test  stand  and  (b)  PSV-400-3D  LDV. 


3.4  Mode  Shape  Extraction 

The  PSV  software  computes  a  F  FT  of  the  time  response  data  for  each  DOF  at  every  valid 
measurement  location.  The  user  can  select  to  animate  the  operational  deflection  shapes  (ODS), 
in  th  e  time  d  omain.  A  n  O  DS  is  r  epresentative  o  f  th  e  a  ctual  d  isplacement  o  f  th  e  p  and  at  a 
specified  f requency.  Therefore,  a  11  th  e  mo  des  are  p  resent  in  s  ome  f orm  i  n  an  ODS  w  ith  th  e 
exception  o  f  m  odes  t  hat  ar  e  n  ot  excited.  The  ODS  i  s  de  pendent  of  t  he  i  nput  1  ocation  and 
magnitude  of  the  force  applied.  Any  changes  to  either  input  location  or  magnitude  the  ODS  will 
likely  change.  F  or  example,  if  the  test  article  is  unconstrained  and  it  is  struck  with  too  much 
force  n ear  t he  outer  ed ge  o f  t he  p and,  a  r igid  bod y  m ode  w ill  a ppear  d ominant.  E  ach  ODS 
contains  many  modes  but  is  typically  dominated  by  a  single  mode,  when  the  user  is  selecting  a 
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frequency  near  resonance.  T  he  mode  shapes  are  extracted  from  the  measured  ODSs  using  the 
modal  extraction  software  ME-Scope  VES  5.0. 

ME-scope  imports  FRF  data  directly  from  the  PSV  software  through  a  proprietary  Vibrant 
Technologies  file  format.  ME-Scope’s  software  utilizes  curve  fitting  algorithms  to  generate  an 
FRF.  A  t  ransfer  f  unction  i  s  c  alculated  b  y  a  ssigning  pol  es  a  nd  z  eros  t  o  t  he  F  RF  c  omer 
frequencies.  The  calculated  transfer  function  fits  a  similar  curve  to  the  measured  FRF.  Natural 
frequencies,  da  mping,  and  m  ode  s  hapes  a  re  t  hen  e  stimated  f  rom  t  he  t  ransfer  f  unction.  T  he 
modes  are  independent  of  the  loads  imparted  on  the  structure  during  the  test.  ME-scope  can  also 
animate  the  modes  by  scaling  the  eigenvectors  with  a  sinusoidal  magnitude  function  in  the  time 
domain.  Animating  these  different  modes  allows  for  a  visualization  of  each  particular  node  and 
erroneous  data  can  easily  be  discovered  and  eliminated.  T ypically,  these  bad  data  points  arrive 
from  scattering  of  the  laser  light  resulting  in  a  poor  return  to  the  LDV  sensor.  A  s  mentioned 
earlier,  a  light  scattering  material  which  is  typically  used  when  scanning  metallic  surfaces  wasn’t 
used  because  PnPSAT  is  flight  hardware.  This  led  to  the  laser  hitting  scratches  in  the  metallic 
surface  leading  to  some  erroneous  data;  needless  to  say,  all  3  D  OF  for  each  bad  data  point  was 
dismissed.  Due  to  the  scattering  effect,  30%  and  18%  of  the  data  points  for  the  -Zand  -Y 
panels,  r  espectively,  were  r  emoved.  After  eliminating  ha  d  da  ta  poi  nts  a  nd  e  xtracting  t  he 
complex-valued  modes,  the  next  step  is  to  analyze  the  data. 

Before  t  he  ex  tracted  m  ode  s  hapes  can  b  e  u  tilized  for  F  E  m  odel  t  uning,  t  he  extracted 
complex-valued  m  ode  s  hape  m  atrix  0^  is  ap  proximated  t  o  a  r  eal-valued  mode  m  atrix  0^;  by 


Neidbal’s  [1984]  complex  transformation: 


where  all  values  are  defined  in  Nomenclature.  Based  on  pr  evious  research  efforts,  Doupe  et  al 
[2009],  it  was  found  Neidbal’s  complex  transformation  yielded  the  best  results  for  this  case. 

3.5  Tuning -Y  and  -Z panels 


After  converting  the  complex-valued  modes  to  real-valued,  the  tuning  process  can  begin.  The 
optimization  process  will  adjust  various  stiffness  related  parameters  in  the  FE  model  until  the 
natural  frequencies  are  tuned  to  be  within  a  predetermined  tolerance  and  the  difference  between 
the  m  ensured  and  analytical  m  odes  i  s  m  inimized.  Again,  t  he  obj  ective  f  unction  that  is 
minimized  is. 


P  2  P  ''  S 


at  each  optimization  step  the  eigenvalues  and  modes  are  compared  to  the  measured.  For  tuning 
the  panels  a.  =  0  and  b.  =  I ;  meaning  the  objective  function  minimizes  the  difference  between 


the  modes  while  constraining  the  natural  frequencies  separately  with  a  t  olerance  given  by  the 
analyst.  Nastran  uses  a  g  radient-based  optimization  method  in  order  to  minimize  the  obj  ective 
function.  A  gradient-based  approach  calculates  partial  derivatives  of  the  objective  function  and 
minimizes  the  objective  f  unction  in  a  f  ew  s  teps  a  s  pos  sible  b  y  adjusting  t  he  pr  eviously 
mentioned  model  parameters  in  the  direction  of  steepest  descent. 

The  first  step  is  to  generate  the  mass  matrix  which  remains  constant  throughout  the  tuning 
process.  In  order  to  calculate  the  reduced  mass  matrix  a  G  uyan  reduction  is  completed.  T  he 
Guyan  reduction  simply  takes  the  measured  DOFs,  which  are  the  basis  for  the  tuning  process, 
and  extracts  them  from  the  FE  model.  The  Guyan  reduction  is  described  by: 
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where  K  is  the  stiffness,  Xj  are  the  measured  DOF  and  are  all  other  DOF.  F  is  the  applied 
load  and  for  this  case  it  is  zero.  From  this,  solve  the  linear  system  of  equations  to  arrive  at 

X2  =  .  (22) 


The  Guyan  reduction  matrix  is: 


G  = 


^22^21 


(23) 


where  I  is  the  identity  matrix.  The  Guyan  reduction  matrix  can  be  applied  to  the  mass  matrix. 

M^=G^MG  (24) 

This  operation  will  reduce  the  mass  matrix  to  the  size  of  the  measured  DOFs.  T  his  matrix  is 
important  for  the  cross-orthogonality  check. 

In  o  rder  toe  nsure  c  onvergence  a  nd  a  r  easonable  s  olution,  c  onstraints  a  re  put  on  t  he 
parameter.  T  ypically  Y  oung’s  modulii  are  allowed  t  o  v  ary  up  to  +/-  20%  of  the  initial  table 
values.  The  natural  frequencies  have  a  fixed  constraint  of  +/-  5%  from  the  measured. 

Mode  t  racking  i  s  a  Iso  us  ed  i  n  or  der  toe  nsure  t  hat  t  he  m  odes  remain  i  n  th  e  in  itial  o  rder. 
Mode  t  racking  i  s  don  e  w  ith  a  cr  oss-orthogonality  ch  eck  b  etween  the  m  ass  n  ormalized 
eigenvectors  from  the  current  and  previous  optimization  steps  and  the  mass  matrix,  represented 
by  the  equation  below: 

(25) 
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where  /  is  the  current  design  cycle  and  t.  is  a  «  by  n  matrix.  If  the  eigenvectors  do  not  change 


between  o  ptimization  c  ycles,  th  e  t,.  will  be  th  e  identity  ma  trix.  If  t  he  eigenvectors  change 

between  optimization  cycles,  some  of  the  values  on  t  he  diagonal  will  be  less  than  1.0.  If  the 
mode  order  has  changed  indicated  by  any  diagonal  value  falling  below  0.9  Nastran  then  corrects 
the  mode  order  even  through  the  associated  natural  frequencies  will  be  out  of  order.  This  mode 
tracking  function  is  important  because  the  optimizer  must  correctly  match  eigenvectors  with  the 
target  modes. 

3.6  Creating  Entire  PnPSAT I  Using  -Y  and  -Z panel 

After  tuning  t  he  -Y  an  d  -Z  p  anels,  the  co  mplete  P  nPSAT  I  m  odel  can  b  e  as  sembled.  In 
FEMAP,  both  the  -Y  and  -Z  panels  are  imported  into  a  single  FE  model  and  joined  together.  The 
panels  are  connected  with  blocks;  a  view  of  two  of  PnPSAT  I  panels  open  and  connected  is  in 
Figure  22. 


Figure  22:  PnPSAT  I  panels  open  and  joined  together  via  hinges. 
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Because  the  +Y,  +X,  -X  panels  are  identical  to  the  -Y  panel,  the  tuned  -Y  panel  FE  model 
can  be  copied  to  create  -Y,  +X,  -X  panel  FE  models.  The  +Z  is  identical  the  -Z  panel,  the  tuned 
-Z  panel  FE  model  can  be  copied  to  create  the  -Z  panel  FE  model.  Two  components  influential 
on  the  tuning  process  of  PnPSAT  I  were  the  base  plate  and  SPA  adapter  ring.  The  base  plate  was 
acquired  from  the  high-fidelity  Spaceworks  model.  Figure  23(b),  and  the  adapter  ring.  Figure 
23(a),  was  created  from  3  adjacent  plate  elements.  The  bottom  of  the  adapter  ring  is  constrained 
in  all  degrees  of  freedom. 


Figure  23:  (a)  Adapter  ring  modeled  with  plate  elements,  constraint  forces  attached  (b)  base  plate 
acquired  from  Spaceworks  model. 


Initially,  the  panels  were  to  be  connected  using  spring  elements;  however,  when  the  panels 
were  joined  together  the  nodes  on  the  side  panels  became  coincident;  meaning  nodes  from  each 
panel  became  a  single  node  sharing  elements  from  both  panels.  A  coincident  node  check  of  the 
entire  s  atellite  m  erged  t  he  n  odes  f  undamentally  r  igidly  attaching  t  he  p  anels  t  ogether.  T  he 
coincident  nodes  can  be  visualized  in  Figure  24(a).  The  SPA  adapter  ring  has  constraint  forces 
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on  it.  These  constraint  forces  constrain  the  three  translational  and  three  rotational,  making  each 
node  constrained  in  all  six  DOFs.  T  he  constrained  adapter  ring  is  attached  via  spring  damper 
elements  to  the  base  plate  (Figure  24(c)).  T  hese  spring  damper  elements’  stiffness  parameters 
will  be  t  uned  toe  nsure  t  he  F  E  m  odd  i  s  a  s  c  lose  t  o  t  he  experimental  n  atural  frequencies  as 
possible. 


Figure  24:  (a)  Demonstration  of  coincident  nodes  (b)  constructed  PnPSAT  I  (c)  side  view  of  the 
spring  elements  connecting  the  adapter  ring  to  the  base  plate. 
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The  b  ase  p  late  is  rigidly  attached  to  th  e  b  ottom  o  f  -Z  p  anel.  T  he  c  onstructed  s  atellite 
structure  can  be  seen  in  Figure  24(b).  The  components  of  the  satellite  to  include  antennas,  GPS 
receivers,  solar  panels,  batteries  are  modeled  as  single  point  mass  elements.  The  components  are 
defined  b  y  t  heir  m  ass  and  the  p  rinciple  M  01  w  hich  cam  e  d  irectly  from  t  he  S  paceworks  FE 
model.  A  quick  check  to  make  the  sure  the  components  are  the  correct  mass  and  in  the  right 
place  is  check  total  mass  and  center  of  gravity  (C.G.).  E  ach  element  is  rigidly  linked  to  nodes 
near  it,  in  essence  providing  an  accurate  yet  more  simplified  model.  A  look  into  the  satellite  at 
the  components  on  the  interior  is  provided  in  Figure  25. 


Figure  25:  Look  into  PnPSAT  I  FE  model  showing  interior  components  attached  via  rigid  links. 


The  complete  PnPSAT  I  FE  model  is  in  Figure  26. 
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Figure  26:  Complete  PnPSAT  I  model  with  components  connected  by  rigid  links. 


3, 7  Measurement  Method  for  Frequency  data  for  Entire  PnPSAT  I 

Now  that  both  the  panels  are  correctly  tuned  and  the  untuned  full  scale  model  is  created  the 
next  step  is  to  tune  the  entire  satellite.  The  method  used  to  conduct  the  forced  response  testing  or 
vibration  of  the  entire  satellite  was  conducted  at  Kirtland  AFB  in  Albuquerque,  New  Mexico.  A 
shaker  table,  as  shown  in  Figure  27,  applied  a  burst  random  vibration  and  burst  sine  ranging  from 
0  to  2000  Hz.  This  shaker  table  is  considered  electro  hydraulic;  with  hydraulic  shakers,  force  is 
generated  through  the  use  of  hydraulics,  which  can  provide  much  higher  force  levels  than  other 
shakers.  The  best  choice  of  excitation  function  depends  on  several  factors  such  as  the  available 
signal  processing  equipment,  characteristics  of  the  structure,  general  measurement  considerations 
and,  of  c  ourse,  t  he  e  xcitation  s  ystem.  In  t  his  c  ase,  bur  st  r  andom  vi  bration  a  nd  bur  st  s  ine 
vibration  w  ere  c  hosen.  H  owever,  ot  her  va  rious  e  xcitation  f  unctions  a  nd  m  ethods  e  xist  a  nd 
possibly  better  exploited  given  the  differences  in  factors  mentioned  above. 
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Now  that  the  methods  for  testing  have  been  described,  the  method  for  acquiring  the  data  for 
this  experiment  must  be  explained.  In  this  case,  accelerometers  were  placed  on  the  satellite  in 
several  different  places.  These  are  noticeable  in  Figure  27  represented  by  the  wires  coming  off 
the  satellite  hanging  down  from  the  shaker  table. 


Figure  27:  PnPSAT  I  on  the  vibration  table  on  Kirtland  AFB. 

The  satellite  is  tested  in  all  three  axes.  E  ach  accelerometer  is  aligned  with  each  axis.  W ith 
this  data  and  using  the  process  described  in  Section  2.4,  the  FRFs  of  each  different  axis  case  can 
be  computed.  Examining  the  FRF  it  is  easy  to  estimate  natural  frequencies.  The  first  few  natural 
frequencies  from  t  he  t  hree  d  ifferent  tests  are  u  sed  for  t  uning.  H  owever,  t  he  m  ode  s  hapes  of 
PnPSAT  I  were  not  taken  and  therefore  the  untuned  model  will  necessarily  be  representing  the 
correct  mode  shapes  as  the  actual  satellite  would  behave.  An  example  of  the  accelerometer  data 
taken  during  the  test  of  the  X  axis  test  is  given  in  Figure  28;  also,  a  comparison  between  the  X 
axis  random  vibration  sweep  and  sine  sweep  FRF. 
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(b) 

Figure  28:  (a)  FRF  of  Sine  sweep  from  0  to  2000  Hz  on  X  Axis  test  (b)  FRF  of  random  vibration 
sweep  from  0  to  2000  Hz  on  X  Axis  Test,  green  represents  input  and  black  represents  output. 
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3,8  Tuning  PnPSAT I 


PnPSAT  I ’s  tuning  process  involved  minimizing  difference  b  etween  the  first  three  KAFB 
measured  natural  frequeneies  and  the  FE  model  natural  frequencies  to  within  +/-  5%.  The  modal 
data  was  not  extracted  from  the  KAFB  data  due  to  the  limited  number  of  aceelerometers,  so 
mode  traeking  is  not  required;  because  the  panel  stiffness  and  mass  properties  were  previously 
tuned  as  explained  earlier  in  this  document.  The  -X,  +X,  -Y,  +X,  -Z,  Z  panel’s  material  eards 
remain  constant  throughout  the  tuning  proeess.  The  6  D OF  spring/damper  element’s  stiffness 
and  the  base  plate  and  adapter  ring’s  Young’s  modulii  serve  as  the  tuning  parameters.  There  are 
a  total  of  148  spring/damper  elements  in  Eq.  ( 19)  in  the  model  that  can  be  tuned,  where  a.  =  1 

and  b.=0.  A  ttaehed  are  18  property  cards;  each  property  eard  possesses  6  D  OF  allowing  for 

1 08  de  sign  va  riables.  An  a  dditional  20  de  sign  va  riables  a  re  t  he  a  dapter  r  ing  a  nd  ba  se  pi  ate 
Young’s  modulii,  totaling  128  design  variables.  The  design  variables  are  altered  by  Nastran  in  a 
gradient-based  optimization  approach  in  order  to  minimize  the  objective  funetion.  Convergence 
is  a  chieved  w  hen  t  he  o  bjective  f  unction  c  annot  be  m  inimized  f  urther  by  manipulation  of  t  he 
design  v  ariables.  Constraints  a  re  a  pplied  on  t  he  6  D  OF  s  pring  e  oefficients  b  y  allowing  t  he 
spring/damper  elements  to  be  altered  by  .01%  to  1000%  of  their  original  values.  The  Young’s 
modulii  are  constrained  to  be  within  +/-40%  of  their  original  table  values.  In  the  next  chapter  we 
discuss  the  results. 
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IV.  Results  and  Discussion 


In  this  chapter  the  tuning  results  of  the  -Y  panel  and  -Z  panel,  and  satellite  FE  model  are 
presented.  In  each  section,  the  untuned  FE  models  are  presented  followed  by  a  comparison  of 
the  tuned  model  to  the  experimental  data. 


4.1  Untuned  —Y  and  -Z panel  FE  model  Results 


First,  we  look  at  the  untuned  panel  FE  model  modes  and  natural  frequencies  to  ensure  the 
experimental  an  d  F  E  m  odels  m  odes  a  gree  w  ith  respect  t  o  each  o  ther.  In  or  der  t  o  a  ssume  t  he 
mass  matrix  is  correct  the  mass  of  the  FE  model  must  match  to  within  a  1%.  It  is  important  that 
the  FE  model  is  in  close  agreement  with  measured  mass.  The  difference  between  the  measured 
and  m  odel  m  ass  s  hould  be  w  ithin  1%  .  A  comparison  be  tween  m  ensured  a  nd  a  nalytical  i  s 
summarized  in  Table  1. 


Table  1 :  Comparison  between  actual  and  model  mass. 


Measured  (kg) 

Analytical  (kg) 

%  Difference 

-Z  Panel 

6.72 

6.7077 

0.183035714 

-Y  Panel 

6.82 

6.8589 

0.570381232 

The  untuned  -Y  and  -Z  panels  have  fewer  FE  DOF  than  the  Spaceworks  FE  panel  models. 
The  comparison  between  the  number  of  nodes  and  elements  for  each  panel  is  summarized  below. 

Table  2:  Comparison  between  Spaceworks  and  newly  created  untuned  -Y  and  -Z  FE  models. 


Spaceworks 
-Z  Panel 

Spaceworks 
-Y  Panel 

Reduced 
-Z  Panel 

Reduced 
-Y  panel 

Nodes 

-42,000 

-48,000 

7,647 

8,875 

Elements 

-28,000 

-30,000 

4,445 

5,273 

It  is  important  to  ensure  no  mechanisms  or  localized  modes  exist  in  the  untuned  FE  models. 
An  eigenvalue  and  normal  mode  shape  test  is  run  on  t  he  untuned  models.  T  his  test  gives  the 
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natural  frequencies  along  with  t  he  c  orresponding  m  odes.  Looking  t  hrough  t  he  m  ode  s  hapes 
calculated  it  is  important  to  ensure  that  the  experimental  and  analytical  modes  agree  in  order  for 
the  t  uning  pr  ocess.  Upon  i  nspection  of  t  hese  m  odes  it  w  as  e  vident  th  at  lo  cal  mo  des  e  xisted. 
These  local  modes  were  caused  by  elements  not  being  properly  linked  to  each  other,  typically  the 
rigid  link  elements  which  connect  the  hinges  on  each  panel.  T  o  fix  this  problem  all  the  rigid 
links  in  the  model  were  erased.  The  second  time  around  it  was  critically  important  to  make  the 
hinge  node  be  dependent,  while  making  the  panel  itself  the  independent  nodes.  This  ensured  that 
whatever  the  nodes  on  t  he  panel  were  doing  the  dependent  nodes  would  mirror.  S  uccessfully 
eliminating  t  his  1  ocalized  m  ode  pr  oblem  t  hereby  allowing  f  or  t  he  m  odes  t  o  1  ine  up.  T  he  s  ix 
mode  shapes  for  both  panels  will  be  compared  to  the  six  measured  mode  shapes.  Refer  to  Figure 
29  for  t  he  -Y  p  anel  an  d  Figure  30  for  t  he  -Z  P  anel;  note,  t  he  m  easured  m  ode  s  hapes  are 
interpolated  from  the  points  measured  by  the  scanning  LDV. 
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Measured 


Mode  1 

235  Hz  Measured 
258  Hz  Analytical 


Mode  2 

281  Hz  Measured 
347  Hz  Analytical 


Mode  3 

346  Hz  Measured 
402  Hz  Analytical 


Mode  4 

461  Hz  Measured 
549  Hz  Analytical 


igi 


Mode  5 

585  Hz  Measured 
662  Hz  Analytical 


Mode  6 

612  Hz  Measured 
689  Hz  Analytical 


Figure  29:  Measured  versus  analytieal  mode  shapes  -Y  panel. 


54 


Measured 


Analytieal 


Mode  1 

26 1  Hz  Measured 
256  Hz  Analytical 


Mode  2 

296  Hz  Measured 
340  Hz  Analytical 


Mode  3 

390  Hz  Measured 
506  Hz  Analytical 


Mode  4 

555  Hz  Measured 
614  Hz  Analytical 


Mode  5 

681  Hz  Measured 
694  Hz  Analytical 


Mode  5 

538  Hz  Measured 
738  Hz  Analytical 


Figure  30:  Measured  versus  analytieal  mode  shapes  -Z  panel. 
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4.2  Tuned -Y  and  -Z panel  FE  model  Results 


The  tuning  process  had  to  be  re-accomplished  several  times  to  get  acceptable  results.  First, 
there  was  an  insufficient  number  of  tuning  parameters  which  resulted  in  failed  tuning  cases.  It 
was  necessary  to  go  back  into  the  FE  model  and  increase  the  number  of  material  cards  present. 
Table  3  summarizes  the  before  and  after  material  cards 


Table  3:  Material  cards  -Z,  -Y  panels. 


-Z  Panel 

-Y  Panel 

Initial  Material  Cards 

14 

10 

Final  Material  Cards 

65 

95 

The  o  bj  ective  function  a  s  w  ritten  is  v  ery  s  ensitive  to  s  mall  ta  rget  v  alues  b  ecause  th  e  ta  rget 
values  are  in  the  denominator  of  J,  these  small  values  make  the  obj  ective  function  erroneously 
large.  This  makes  the  tuning  process  very  sensitive  to  only  a  few  DOF  in  each  mode.  After 
discarding  these  data  points  with  very  small  values  and  plotting  the  normalized  untuned  modes 
by  t  he  m  ensured  i  t  b  ecame  cl  ear  t  he  o  bj  ective  f  unction  w  as  ab  solved  oft  hese  f  ew  bad  d  ata 
points.  A  n  example  of  a  bad  obj  ective  function  and  a  correct  obj  ective  function  for  the  first 
mode  of -Z  panel  can  be  seen  in  Figure  3 1 . 


Analytical 
Divided  by 
Measured 


(a)  (b) 

Figure  31:  -Z  Panel  (a)  Unsatisfactory  objective  function,  single  value  approaching  300,  tuned 
mode  di  vided  b  y  s  mall  e  rroneous  da  ta  ( b)  R  emoved  ba  d  d  ata  poi  nts  a  11  poi  nts  w  ith  +  /-  3 
magnitude  -  note  the  y-axis  scale  difference. 
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With  all  the  bad  nodes  removed,  a  visual  check  shows  that  all  the  nodes  are  in  their  correct 
locations  and  plane.  There  were  130  and  137  measurement  locations  in  the  -Z  and  -Y  panels, 
respectively.  The  nodes  being  tuned  are  represented  by  circles  in  Figure  32: 


Figure  32:  (a)  -Z  panel,  130  nodes  (b)  -Y  panel,  137  nodes. 
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Tuning  results  for  the  -Y  and  -Z  panel  are  loeated  in  Table  4  and  Table  5.. 


Table  4:  Measured  versus  tuned  for  modes  1-6  for  -Y  panel. 


Measured  vs  Untuned 

Measured  vs  Tuned 

Mode 

Measured  (Hz) 

Untuned  (Hz) 

Abs%DifF 

Tuned  (Hz) 

Abs  %  Diff 

1 

235.4 

258.0 

9.60 

229.5 

2.51 

2 

281.6 

347.0 

23.22 

289 

2.63 

3 

346.1 

402.0 

16.15 

348.8 

0.78 

4 

461.3 

549.0 

19.01 

472.8 

2.49 

5 

585.2 

662.0 

13.12 

574.9 

1.76 

6 

612.2 

689.0 

12.54 

597.6 

2.38 

Table  5:  Measured  versus  tuned  for  modes  1-6  for  -Z  panel;  ordered  by  analytical  modes. 


Measured  vs  Untuned 

Measured  vs  Tuned 

Mode 

Measured  (Hz) 

Untuned  (Hz) 

Abs  %  Diff 

Tuned  (Hz) 

Abs  %  Diff 

1 

261.2 

256.0 

1.99 

249.2 

4.59 

2 

296 

340.0 

14.86 

294.3 

0.57 

3 

390 

506.0 

29.74 

401.2 

2.87 

4 

555.5 

614.0 

10.53 

528.4 

4.88 

5 

680.7 

694.0 

1.95 

564.3 

17.10 

6 

538 

738.0 

37.17 

648.8 

20.59 

Observing  Table  4  it  is  evident  that  -Y  panel  was  tuned  to  be  within  +!-  3%  of  the  measured 
data;  however,  the  -Z  panel  proved  much  more  difficult  to  tune.  The  FE  model  had  modes  in  a 
different  order  than  the  measured  data  had  them  being  shown.  The  5*  and  6*  measured  natural 
frequencies  are  6 80.7  and  5  38  Hz,  respectively.  The  measured  and  analytical  modes  agree  in 
terms  of  order,  but  the  tuning  process  had  poor  agreement  with  the  frequencies.  After  comparing 
the  f  requencies  i  t  w  as  i  mportant  to  1  ook  at  t  he  comparison  b  etween  an  alytical  and  m  easured 
modes.  U  sing  a  cross  orthogonality  check  between  the  mass  normalized  measured  and  tuned 
eigenvectors  we  arrive  at  a  value  from  0  to  1,  1  b  eing  nominal.  A  comparison  of  the  measured 
versus  measured,  and  measured  versus  tuned  are  in  Figure  33  with  a  comparison  of  the  cross- 
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orthogonality  terms  in  Table  6.  The  comparison  is  done  using  a  mass  weighted  orthogonality 


check.  Given  as  Eq.  (26)  were  is  the  Guyan  reduced  mass  matrix. 


(26) 


Experimental  Modes 


(a) 


Analytical  Modes 


Experimental  Modes 


Experimental  Modes 


(b) 


Figure  33:  (a)  -Y  panel  tuned  cross-orthogonality  results  (b)  -Y  measured  vs.  measured  cross¬ 
orthogonality  (c)  -Z  panel  tuned  cross-orthogonality  results  (d)  -Z  measured  vs.  measured  cross¬ 
orthogonality 
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Table  6:  Table  comparing  untuned  and  tuned  cross  orthogonality  checks  -Z  panel. 


Mode  # 

XO  Untuned 

XO  Tuned 

1 

0.995 

0.996 

2 

0.970 

0.976 

3 

0.935 

0.949 

4 

0.200 

0.204 

5 

0.937 

0.941 

6 

0.895 

0.865 

Table  7:  Table  comparing  untuned  and  tuned  cross  orthogonality  checks  -Y  panel. 


Mode  # 

XO  Untuned 

XO  Tuned 

1 

0.977 

0.973 

2 

0.947 

0.229 

3 

0.929 

0.923 

4 

0.895 

0.870 

5 

0.102 

0.952 

6 

0.879 

0.822 

Observing  the  measured  versus  measured  cross  orthogonality  plots  it  is  important  to  notice 
the  similarities  between  the  4th  and  5th  mode  for  the  -Z  panel  and  the  2nd  and  5th  mode  for  the  - 
Y  panel.  For  each  panel  the  similarities  that  exist  are  noticeable  on  the  cross  orthogonality  plots. 
These  modes  being  so  closely  aligned  affect  the  tuning  ability  because  the  optimization  process 
has  difficulty  determining  the  difference  between  these  modes. 

A  brief  look  at  the  tuned  Young’s  modulii  for  the  -Y  and  -Z  panel  are  shown  in  Figure  34 
and  Figure  35,  respectively.  The  initial  value  of  the  Young’s  modulii  for  both  models  was  9.9E6 
Ksi. 
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(a) 


Figure  34:  (a)  Modulus  of  Elasticity  of  the  various  tuning  variables  the  -Y  Panel.  Initial  value  of 
9.9E6  Ksi.  (b)  E  lements  w  ith  5.90E  6  modulus  of  e  lasticity.  (  c)  E  lements  w  ith  m  odulus  of 
elasticity  b  etween  5.90E  6  K  si  a  nd  1 .30E  7  K  si.  ( d)  E  lements  w  ith  m  odulus  of  e  lasticity  of 
1.30E7  Ksi. 
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Figure  35:  (a)  Modulus  of  Elasticity  of  the  various  tuning  variables  the  -Z  Panel.  Initial  value  of 
9.9E6  K  si.  ( b)  E  lements  w  ith  3.50E  6  modulus  of  e  lasticity.  (  c)  E  lements  w  ith  m  odulus  of 
elasticity  b  etween  3.50E  6  K  si  a  nd  1 .60E  7  K  si.  ( d)  E  lements  w  ith  mo  dulus  o  f  e  lasticity  o  f 
1.60E7  Ksi. 
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4.3  Untuned PnPSAT I 


The  -Y  and  -Z  panels  deseribed  were  ehosen  to  represent  the  -X,  +X,  Y  and  Z,  respeetively, 
for  t  he  s  tructure  of  P  nPSAT  I.  T  he  c  ompleted  s  tructure  with  c  omponents  F  E  m  odel  w  as 
approximately  six  times  smaller  in  size  than  the  Spaeeworks  model.  The  differenee  in  size  can 
be  seen  in  Table  8. 

Table  8:  Comparison  of  size  of  Spaeeworks  and  reduced  PnPSAT  I FE  models. 


Spaeeworks 
PnPSAT  I 

Reduced 

PnPSAT  I 

Factor 

Smaller 

Nodes 

368,850 

67,381 

5.5 

Elements 

246,295 

40,861 

6.0 

DOF 

1,104,633 

208,404 

5.3 

With  PnPSAT  I  fully  constructed  the  mass  is  validated  to  be  within  a  3%  tolerance  window. 
The  a  ctual  m  ass  of  P  nPSAT  I  w  as  1 13.143k  g,  t  he  FE  m  odel  w  as  1 1  5.854,  t  his  i  s  a  2.4  % 
difference;  w  ell  within  th  e  to  lerance  limit.  An  e  igenvalue  and  eigenvector  analysis  was 
conducted  to  ensure  no  localized  modes  or  singularities  exist  before  the  tuning  process  began. 
The  first  natural  frequency  and  mode  shape  analysis  revealed  the  natural  frequencies  were  to  low 
and  the  base  plate  and  -Z  panels  were  not  adequately  fastened  with  rigid  links.  Increasing  the 
natural  frequency  was  accomplished  by  conducting  a  coincident  node  check  to  rigidly  attach  the 
panels  together.  To  fasten  the  base  plate  to  the  -Z  panel,  more  rigid  links  were  added. 

The  parameters  to  be  altered  during  the  optimization  process  are  the  spring/dampers  which 
attach  the  base  plate  and  the  adapter  ring.  This  is  justified  because  the  panels  have  already  been 
tuned  and  the  base  plate  and  adapter  ring  are  the  main  sources  of  strain  on  the  satellite.  Looking 
at  the  strain  energy  for  the  first  mode  in  Figure  36  it  is  evident  that  the  Young’s  modulii  and 
stiffness  o  f  th  e  s  pring/dampers  is  mo  st  imp  ortant  d  uring  th  e  c  omplete  PnPSAT  I  FE  mo  del 
tuning  process. 
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Figure  36:  Blown  up  vi  ew  be  tween  b  ase  pi  ate  and  a  dapter  r  ing,  p  rimary  strain  a  long  t  op  of 
adapter  ring. 


The  initial  stiffness  values  are  allowed  to  change  between  a  tolerance  of  .01%  and  1000%  of 
the  starting  value.  T  he  Y oung’s  modulii  for  the  base  plate  and  adapter  ring  are  allowed  to  be 
altered  by  up  to  +/-40%  of  the  initial  value.  The  first  three  modes  and  natural  frequencies  for  the 
untuned  PnPSAT  I  are  seen  in  Figure  37,  Figure  38,  Figure  39. 
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Figure  37:  Mode  1:  Rocking  X:  untuned  63.2  Hz:  (a)  undeformed  dimetric  view  (b),  (c)  scaled 
up  deflections  demonstrating  the  rocking  motion. 


Figure  38:  Mode  2:  Rocking  Y:  untuned  70.6  Hz  (a)  undeformed  dimetric  view  (b),  (c)  scaled 
up  deflections  demonstrating  the  rocking  Y  mode. 
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4.4  Tuned PnPSAT I 


The  tuning  process  for  PnPSAT  I  was  conducted  with  ease  after  the  experience  gained  from 
tuning  the  -Y  and  -Z  panels.  The  optimization  process  using  4  dual  cores  AMD  8200-2.8  GHz 
with  64  gigabytes  of  RAM  took  approximately  6  hours.  The  tuning  process  length  grows  based 
on  operations,  where  «  is  the  number  of  DOFs  for  the  model.  T he  S paceworks  m odd  i s 
approximately  6  times  larger  and  the  operations  are  cubed  therefore  making  the  total  time  per  run 
unrealistically  long. 

The  tuning  process  for  PnPSAT  I  involved  getting  the  natural  frequencies  for  the  first  three 
modes  to  be  within  5  %  of  experimental  data.  A  table  summarizing  the  difference  between  the 
post  tuned  analytical  and  experimental  is  in  Table  9. 


Table  9:  Experimental  vs.  tuned  natural  frequencies  for  modes  1-3. 


Mode 

Measured  (Hz) 

Untuned  (Hz) 

Abs  %  Diff 

Tuned  PnPSAT  I 
(Hz) 

Abs  %  Diff 

1 

65.00 

63.2 

2.77 

63.8 

1.85 

2 

69.00 

78.5 

13.77 

70.0 

1.45 

3 

190.00 

195.8 

3.05 

194.9 

2.58 

The  table  clearly  shows  the  tuning  process  tuned  the  first  three  natural  frequencies  to  within 
+/-3%  of  the  measured.  The  tuned  Young’s  modulii  values  for  the  tuning  variables  are  in  Figure 
40  while  the  spring  stiffness  coefficients  are  in  Figure  41.  Notice,  the  ceiling  and  floor  were 
reached  by  several  variables  especially  while  tuning  the  springs. 
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5.90E6  5.90  E6  <  E  <  1 .30E7  1 .30E7 

Young’s  Modulus  (Ksi) 


(C)  (d) 

Figure  40:  (a)  Modulus  of  Elasticity  of  the  various  tuning  variables  the  -Y  Panel.  Initial  value  of 
9.9E6  K  si.  ( b)  E  lements  w  ith  5.90E  6  modulus  of  e  lasticity.  ( c)  E  lements  w  ith  m  odulus  of 
elasticity  b  etween  5.90E  6  K  si  a  nd  1 .30E  7  K  si.  ( d)  Elements  w  ith  mo  dulus  o  f  e  lasticity  o  f 
1.30E7  Ksi. 


69 


Figure  41:  Spring  Stiffness  Coefficients:  Includes  both  axial  and  torsional  spring  values.  The 
initial  starting  parameter  value  was  1.0E3  Ibf/in  (axial)  or  1.0E3  Ibf  in/rad  (torisonal). 
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V.  Conclusion 


5.1  Research  Overview 

The  P  nP  c  oncept  on  which  P  nPSAT  I  w  as  i  nitially  conceived  is  aimed  at  t  he 
elimination  or  minimization  of  expensive  and  time  consuming  AIT.  The  ORS  goal  is  to 
launch  payloads  to  space  in  only  six  days.  The  objective  of  this  research  is  to  investigate 
the  p  ossibility  o  f  c  reating  an  accurate  FE  m  odel  o  f  P  nPSAT  I  which  can  accurately 
predict  the  dynamic  response  of  PnPSAT  I. 

In  Chapter  II,  we  discussed  how  to  accomplish  the  process  of  creating  an  FE  model, 
the  methods  used  to  extract  frequency  information  from  an  object  and  the  mathematical 
techniques  us  ed  t  o  c  onduct  m  odal  a  nalysis.  Also  d  iscussed  w  as  FE  an  alysis  which 
included,  the  assumptions  implemented  to  create  an  FE  model,  modal  analysis  from  the 
equations  of  motion  (EOM)  through  the  eigenvalue  problem  (EVP)  leading  to  the  modal 
frequency  response  pr  oblem.  Lastly,  t  he  m  ethod  us  ed  t  o  e  xtract  da  ta  f  or  t  uning,  PS V 
scanning  L  DV,  a  nd  the  obj  ective  function  ut  ilized  f  or  t  uning  t  he  F  E  model  w  ith  t  he 
acquired  mode  shapes  and  natural  frequencies. 

In  C  hapter  III,  w  e  di  scussed  how  t  he  original  P  nPSAT  I F  E  m  odel  w  as  cr  eated  b  y 
Spaceworks  and  contained  1,104,633  D  OFs.  The  size  of  this  model  made  it  extremely 
computationally  expensive.  A  simple  eigenvalue  analysis  on  the  model  took  18  minutes 
on  4  dua  1  c  ore  A  MD  8  200-2. 8GHz  64  g  igabytes  of  R  AM  m  achine.  The  e  igenvalue 
problem  is  solved  thousands  of  times  during  the  tuning  process  therefore  the  size  of  the 
PnPSAT  I  FE  model  needed  to  be  reduced  in  order  to  allow  the  tuning  process  to  finish  in 
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a  reasonable  amount  of  time.  The  creation  of  a  reduced  FE  model  was  started  when  the 
ORS  office  agreed  to  allow  ping  testing  on  the  -Y  and  -Z  panels.  Assuming  the  -X,  +X, 
Y  and  Z  panels  are  the  same  as  the  -Y  and  -Z  panels,  respectively,  it  w  as  possible  to 
generate  an  accurate  representation  of  PnPSAT  I  out  of  two  correctly  tuned  panels.  From 
this  the  geometry  of  the  -Y  and  -Z  panels  were  extracted  from  PnPSAT  Fs  CAD  file  the. 
A  ne  twork  of  node  s  was  c  reated  b  y  h  and  i  n  F  EM AP  a  nd  meshed  t  ogether  b  y  ha  nd  t  o 
create  a  fluid  FE  model  for  both  the  -Y  and  -Z  panels. 

Ping  testing  using  the  PSV  scanning  LDV  was  conducted  to  collect  modes  and  natural 
frequencies  for  bot  h  t  he  -Y  an  d  -Z  p  anels.  This  w  as  accomplished  us  ing  a  n  M  B 
Dynamics  C  al50  E  xciter  el  ectrodynamic  s  haker  b  eing  d  riven  b  y  a  H  ewlett  P  ackard 
33120A  15M Hz/Arbitrary  waveform  generator  p rogrammed  t o  g enerate  a  burst  square 
wave  t  o  ex  cite  t  he  p  anels.  T  he  t  ests  p  rovided  e  xcellent  d  ata  i  n  the  1  ower  f  requencies, 
below  ~  500  H  z,  w  ith  excellent  co  herence  and  s  pectral  d  ensities.  T  he  ex  perimentally 
gathered  d  ata  w  as  t  hen  s  pectrally  s  ieved  u  sing  M  E-Scope  V  ES  5  .0  t  o  acq  uire  clean 
averaged  frequency  response  data.  The  tuning  process  involved  ensuring  analytical  and 
experimentally  g athered  mode  differences  are  minimized  while  ensuring  the  analytical 
and  experimentally  collected  natural  frequencies  are  held  within  their  constraint  limits  of 
+/-5%. 

The  t  uned  p  anel  m  odels  ar  e  co  pied  an  d  at  tached  t  o  cr  eate  an  i  nitial  s  tructural  F  E 
model.  The  -Z  panel  and  base  plate  are  rigidly  attached  via  rigid  links.  The  adapter  ring 
and  the  base  plate  are  attached  via  6  DOF  spring/damper  elements.  Lump  masses  with 
correct  ma  ss  a  nd  mo  ment  o  f  in  ertia  in  a  11  th  ree  d  irections  a  re  r  igidly  a  ttached  to  th  e 
structure  leading  to  a  complete  PnPSAT  I  FE  model.  This  complete  FE  model  is  tuned 
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once  b  y  ensuring  t  he  e  xperimentally  gathered  and  an  alytical  n  atural  f  requencies  ar  e 
within  +1-5%  of  one  another  for  the  first  three  modes. 

In  Chapter  IV,  we  discussed  the  -Y  and  -Z  panels  were  6.3  and  5.7  times  smaller  than 
the  Spaceworks  -Y  and  -Z  panels  respectively.  The  -Y  and  -Z  panels  were  tuned  to  be 
within  +1-5%  and  20%  of  the  first  six  experimentally  collected  natural  frequencies  and 
modes,  respectively.  T  he  Y oung’s  modulii  were  allowed  to  fluctuate  up  to  +/-40%  of 
their  original  values.  The  -Y  panel  was  tuned  within  3%  for  all  six  natural  frequencies. 
The  -Y  panel  mode  tuning  are  all  within  20%  with  the  exception  of  the  5*  mode.  The 
experimentally  collected  data  for  the  2”^*  and  5*’’  mode  are  nearly  the  same.  The  tuning 
algorithm  h  as  d  ifficulties  d  isceming  b  etween  th  e  tw  o  th  erefore  d  uring  th  e  c  ross 
orthogonality  checks  the  2”^*  mode  gets  shoved  down  and  the  5*  mode  is  shoved  up.  The 
-Z  pa  nel  w  as  t  uned  w  ithin  5%  f  or  the  f  irst  4  modes,  t  he  1  ast  t  wo  m  odes  bus  ted  t  he 
constraint  because  the  measured  natural  frequencies  for  the  5*  and  6*’’  mode  are  out  of 
order  with  the  FE  model.  T  he  remedy  for  this  was  to  swap  the  natural  frequencies  and 
modes  in  attempt  for  the  tuning  algorithm  to  adjust  the  models  modes  to  the  correct  order. 
The  tuning  algorithm  was  not  able  to  overcome  this  for  the  natural  frequencies;  however, 
the  m  ode  s  hapes  do  m  eet  t  he  20%  c  onstraint  w  ith  t  he  e  xception  of  t  he  4  mode.  I  n 
similar  fashion  as  the  -Y  panel  the  4*  and  5*’’  mode  are  very  similar  to  one  another;  hence 
why  the  4*  mode  is  not  within  the  constraint  limit. 

The  completed  PnPSAT  I  model  had  approximately  6  times  less  DOF  than  the  larger 
Spaceworks  m  odel.  T  he  P  nPSAT  I  w  as  t  uned  t  o  be  w  ithin  +  /-5%  o  f  th  e  f  irst  th  ree 
experimentally  collected  natural  frequencies.  The  Young’s  modulii  of  the  base  plate  and 
the  ad  apter  ring  w  ere  allowed  t  o  fluctuate  u  p  t  o  +/  -40%  of  t  heir  o  riginal  va  lues.  T  he 
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spring/damper  element’s  stiffness  varied  from  .01%  to  1000%  of  their  original  values. 
The  t  uned  P  nPSAT  I’s  FE  m  odel  natural  f  requencies  are  a  11  w  ithin  in  3  %  of  t  he 
experimentally  gathered  data  from  KAFB  tests. 

5.2  Future  Work 

The  ever  changing  and  complex  space  community  needs  to  be  able  to  become  faster 
and  m  ore  a  gile  i  n  m  eeting  d  emands.  P  nPSAT  I  p  rovides  e  xcellent  groundwork  f  or 
studying  t  he  us  e  of  computationally  i  nexpensive  a  nd  a  ccurate  F  E  m  odels  i  n  or  der  t  o 
predict  dynamic  responses  of  various  satellite  configurations.  Numerous  additional  tasks 
could  be  completed  to  further  this  research.  First,  studies  could  be  focused  on  the  effect 
of  m  aterial  c  ards  on  t  he  t  uning  pr  ocess,  t  o  i  nclude  how  m  any  a  re  i  deal,  w  hat  i  s  a  n 
appropriate  pe  rcentage  the  Y  oung’s  m  odulii  c  an  c  hange.  Second,  a  1  ook  i  nto  t  he 
differences  b  etween  pe  rforming  a  coincident  n  ode  c  heck  v  ersus  us  ing  s  pring/damper 
elements  to  attach  the  panels  and  the  relative  effect  of  both  on  the  tuning  process.  Third, 
a  deeper  look  into  a  more  effective  way  of  connecting  the  point  mass  components  to  the 
satellite  s  hould  be  s  tudied.  T  he  us  e  of  rigid  1  inks  i  s  s  omewhat  a  rbitrary  and  a  s  tudy 
between  different  numbers  of  links  and  locations  of  links  on  stiffness  should  be  explored. 
Fourth,  research  into  the  benefits  of  simply  creating  and  tuning  the  entire  PnPSAT  1  FE 
model  i  nitially;  t  hereby,  s  kipping  t  he  t  uning  p  rocess  o  f  i  ndividual  pa  nels.  Lastly, 
experimentally  co  llected  n  atural  f  requencies  an  d  m  odal  d  ata  co  uld  b  e  acquired  f  or  t  he 
first  10m  odes  o  f  t  he  f ully  c  onfigured  P  nPSAT  1 .  T  his  i  nformation  w  ould  dr  astically 
increase  the  confidence  in  the  tuned  PnPSAT  1  FE  model.  Using  this  tuned  model  it 
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would  be  pos  sible  toe  hange  t  he  e  onfiguration  of  P  nPSAT  I  a  nd  t  he  F  E  m  odel  t  o 
determine  if  an  aeeurately  tuned  structure  can  accurately  predict  dynamic  solutions.  This 
would  he  Ip  m  eet  t  he  O  RS  m  ission  b  y  a  llowing  s  pacecraft  e  ngineers  t  o  pos  sibly  a  void 
AIT  process  all  together  and  allow  them  to  have  full  confidence  in  the  FE  model.  This 
would  ultimately  benefit  the  ORS  mission  by  reducing  the  time  of  AIT  reducing  the  time 
it  takes  to  get  spacecraft  on  orbit  to  benefit  the  war  fighter. 
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